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ABSTRACT 


The  object  of  this  work  is  the  experimental  determination  of 
laminar-turbulent  transition  in  shear  layers  separating  adjacent 
compressible  streams.  An  earlier  theoretical  determination  of  trans¬ 
ition,  based  on  a  minimum  Reynolds  number  condition,  is  reviewed  as  a 
possible  asymptotic  transition  limit  when  the  flow  is  either  an  equilibrium 
shear  layer  or  an  equilibrium  wake.  A  theory  for  the  laminar  flow  devel¬ 
opment  in  the  non-equilibrium  regime  downstream  of  the  trailing  edge  is 
also  presented  for  the  first  time,  and  is  found  to  agree  very  well  with 
numerical  methods.  The  transition  and  laminar  flow  theories  are  then 


joined  to  design  an  experiment  aimed  at  producing  non-equilibrium  flows 
in  a  supersonic  wind-tunnel,  and  at  detecting  the  transition  zone  in  such 
a  flow.  Flow  and  transition  measurements  done  in  this  layer  show  that 
the  shear-layer  transition  Reynolds  number,  based  on  the  layer  thickness 
and  the  fast-side  unit  Reynolds  number,  agrees  with  the  corresponding 
number  found  by  the  equilibrium  transition  theory  within  25%.  It  is 
concluded  that  the  asymptotic  transition  theory  utilized  gives  reasonable 


estimates  even  close  to  the  trailing  edge. 


FOREWORD 


This  report  describes  research  conducted  at  Montana  State 
University's  Supersonic  Wind-Tunnel  Laboratory  under  U.S.  Air  Force 
Office  of  Scientific  Research  Grant  No.  80-0267.  The  program 
consisted  of  a  series  of  experiments  done  to  determine  laminar- 
turbulent  transition  behavior  in  a  two-dimensional  shear  layer 
separating  two  adjacent  supersonic  streams,  with  emphasis  in  the 
non-equilibrium  zone  immediately  downstream  of  the  dividing  partition. 
Responsible  for  this  program  was  A.  Demetriades,  Professor  of  Mechanical 
Engineering.  Timothy  Brower,  graduate  student  in  mechanical  engineering, 
performed  most  of  the  design  work  and  carried  out  most  of  the  measure¬ 
ments. 

The  authors  are  indebted  to  MSU's  John  Rompel ,  Gordon  Williamson, 
and  Glen  McCullough  for  expert  assistance  with  the  laboratory  models 
and  equipment,  and  to  Professor  Harry  Townes  for  his  participation  in 
the  nozzle  design.  Paul  J.  Ortwerth  of  USAF/AFWL  and  Captain  Mike 
Francis  of  AFujR  provided  encouragement  throughout  this  work.  Theore¬ 
tical  work  contained  in  the  Appendix  of  this  report  had  its  inception 
during  earlier  work  of  the  senior  author  for  TETRA  Corp.  of  Albuquerque, 
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1.  Introduction 


It  is  well  known  that  the  mixing  of  two  different  but  parallel, 
adjacent  streams  increases  drasticall;  when  the  shear  layer  separating 
them  becomes  turbulent.  In  fluid  lasers  where  rapid  mixing  of  flowing 
chemicals  is  desired,  the  prediction  of  laminar-turbulent  transition 
onset  is  therefore  of  great  importance.  Mathematical  approaches  leading 
to  such  predictions  are  currently  unavailable.  Traditional  analytic 
methods  center  about  hydrodynamic  stability  theory,  which  however  deals 
with  destabilization  of  flow  disturbances  rather  than  with  turbulence 
onset.  Furthermore,  progress  with  stability  theory  can  be  made  only  in 
inverse  proportion  to  the  complexities  bestowed  on  the  flow  by  design 
features  such  as  compressibility,  heat  transfer  etc.  A  third  obstacle  to 
the  utility  of  stability  theory  is  that  the  shear  layer  profile  is  not 
self-similar  in  the  region  of  interest  because  of  distortions  caused  to 
it  by  the  merging  boundary  layers. 

Apparently,  mixing  can  be  accelerated  by  mechanisms  other  than  turbu¬ 
lence.  The  observations  of  Roshko  and  his  co-worker  (Reference  1)  have 
motivated  researchers  at  Caltech  and  elsewhere  to  seek  inhomogeneous  mixing 
models  based  on  vortex  structures  generated  by  Kel vin-Helmholtz  insta¬ 
bilities.  The  role  of  these  instabilities,  however,  is  not  presently 
clear  when  density  differences  are  also  generated  by  compressibility,  as 
is  the  case  of  the  supersonic  mixing  flows  in  laser  cavities.  While 
supersonic  free  shear  layer  experiments  should  allow  for  observation  and 
assessment  of  such  large-scale  mixing  mechanisms,  they  should  still  aim 
at  disposing  of  the  transition  question  as  the  latter  has  been  understood 
in  the  past. 

A  theory  based  on  minimum  conditions  to  be  satisfied  before  transition 
appears  in  the  layer,  first  proposed  by  the  senior  author  of  this  report 
(Reference  2),  has  been  fairly  successful  when  it  was  compared  with  transi¬ 
tion  data  for  "equilibrium"  free  shear  layers,  i.e.  for  layers  with  an 
asymptotic  self-similar  profile  (Section  2,  below).  There  is  no  reason, 
however,  to  expect  this  approach  to  be  valid  near  the  trailing  edge  of  the 
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partition  separating  the  two  streams,  where  the  shear  layer  is  in  "non¬ 
equilibrium"  caused  by  the  presence  of  the  boundary-layer  "wake".  Experi¬ 
ments  were  thus  performed  in  this  non-equilibrium  region  to  provide  a 
link  to  the  asymptotic  theory.  The  description  of  these  measurements  is 
the  purpose  of  this  report. 
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2.  Predictive  Method  for  Transition  in  Shear  Flows 
2.1  General  Description  of  the  Method 


A  method  for  predicting  transition  to  turbulence  in  free  shear  layers 
was  proposed  by  the  senior  author  in  1979  (References  3  and  4).  A  brief 
description  of  the  principles  behind  this  method  is  most  important  in 
comprehending  the  approach,  results  and  conclusions  of  the  present  experi¬ 
ments. 


Briefly,  the  method  consists  of  establishing  a  minimum  condition  for 
the  maintainance  of  turbulence  in  the  shear  layer,  and  involves  the  follow¬ 
ing  steps.  First,  a  parameter  is  sought  whose  value  in  the  turbulence  must 
lie  above  a  certain  threshold  for  turbulence  to  exist;  this  parameter  is 
the  turbulence  Reynolds  number 

<’) 

based  on  the  r.m.s.  longitudinal  velocity  fluctuation  u',  the  integral  scale 
A  and  the  molecular  kinematic  viscosity  v.  Second,  it  is  claimed  that  if 
self-preserving  turbulence  is  to  follow  transition,  the  quantities  u'  and  A 
must  scale  with  characteristic  velocities  and  lengths  found  in  the  turbulent 
zone,  such  as  the  velocity  defect  in  a  wake  or  the  velocity  "jump”  in  a  FSL, 
and  the  width  of  the  turbulent  zone.  Finally,  the  viscosity  is  expressed 
in  terms  of  the  static  temperature,  changing  eq.  (1)  into  one  containing 
only  mean  (time-averaged)  flow  quantities  such  as  the  critical  flow  width, 
edge  properties,  integral  properties  (e.g.  the  momentum  Reynolds  number) 
and  the  appropriate  empirical  constants.  For  any  specific  type  of  flow 
these  constants  are  known  from  experiments,  so  that  a  set  of  initial  and 
edge  conditions  yield  a  flow  width  value  below  which  the  turbulence  cannot 
sustain  itself. 


It  is  seen  from  the  above  that  the  procedure  gives,  at  best,  a  necessary 
condition  for  turbulence  onset;  no  condition  of  sufficiency  exists.  Threshold 
type  transition  criteria  are  not  new;  for  example  Liepmann  proposed,  long 
ago,  an  analogous  criterion  based  on  turbulent  stress  thresholds  (Reference  5) 
The  attraction  of  such  simple  criteria  is  their  simplicity  of  application, 
so  long  as  they  contain  enough  of  the  physics  of  the  transition  process  and/or 
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the  turbulent  state.  The  physical  notion  contained  in  the  present  method 
is  that  viscosity  damps  the  turbulence;  in  this  respect  the  method  is 
closely  related  to  hydrodynamic  stability  theory.  In  fact  the  proposed 
threshold  of  a  minimum  necessary  Re^  is  analogous  to  the  minimum  critical 
Re0  of  stability  theory,  which  is  used  successfully  to  predict  transition 
trends  in  the  presence,  for  example,  of  pressure  gradients,  cooling  or 
mass  blowing  (Reference  6).  One  application  of  the  present  method  which 
comes  to  mind  immediately,  for  instance,  is  the  minimum  transition  distance 
possible  in  the  presence  of  artificial  disturbances.  A  catalog  of  the 
"classic"  flows  such  as  pipes,  jets  etc.  presented  by  the  senior  author  in 
Reference  4  shows  that  the  observed  effect  of  geometry,  compressibility  and 
heat  transfer  on  these  flows  matches  exactly  the  predictions  of  the  present 
method. 


2.2  Application  to  Make  Flows 


This  qualitative  coincidence  in  trends  induced  this  author  to  see  if 
the  condition  of  necessity  was  also  capable  of  matching  quantitatively 
the  transition  observation  for  various  flows.  An  attempt  was  first  made 
with  two-dimensional  wakes.  In  Reference  7  a  search  of  the  experimental 
literature  was  made  for  the  appropriate  constants,  the  chief  finding  being 
that  the  minimum  possible  value  of  Re^  was 


Other  pertinent  observations  were  that  in  the  "far"  wake  A  -  0,2h,  and  that 
u'  scaled  with  the  velocity  defect  multiplied  by  appropriate  compressibility 
(density)  factors.  With  these  inputs  eq.  (2)  gave 


X  -  XT  .  o(y-l)2  Me^l+Q)2  Re 
at  r  h  " 

- 

1  -1 

rK 

k+1/2  _.j 

T  CDn  1 6* 

V 

(3) 


where  a  and  y  are  the  Prandtl  no.  and  specific  heat  ratio  respectively. 

Me  the  edge  Mach  number,  Q  the  non-dimensional  heating  or  cooling  of  the 
wake,  Cph  the  momentum  defect  (drag)  of  the  wake  and  Rew  the  Reynolds 
number  based  on  it,  k  the  temperature-viscosity  exponent  and  ReQ  is  propor¬ 
tional  to  Re^Q.  Comparison  of  eq.  (3)  with  the  available  data  showed 
satisfactory  agreement  (Reference  7).  A  plot  of  eq.  (3)  appears  on  Figure  1. 
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Quantitative  predictions  for  pipe  flows  and  boundary  layers,  although 
unpublished,  had  also  been  made.  The  transition  Reynolds  number  for  pipes 
was  found  to  be  of  order  of  a  few  thousand,  which  agrees  with  observations. 

For  the  boundary  layer  the  transition  Re.  for  a  smooth  flat  plate  was  found 

a 

to  be  of  order  500  at  low  speeds  and  to  increase  with  Me  beyond  Me  =  4,  to 
increase  with  surface  cooling,  and  to  "reverse"  for  certain  combinations  of 
Me  and  temperature  ratio.  These  trends  agree  with  the  concensus  of  observa¬ 
tions  (and  stability  theory  predictions)  of  boundary-layer  transition, 
although  quantitatively  there  is  desagreement  with  a  fraction  of  existing  data. 


)lication  to  Free  Shear  Layers 


Quantitative  estimates  of  FSL  transition  were  made  in  References  2,  3 
and  4,  following  the  same  procedure  which  led  to  eq.  (3).  As  before,  it 
was  assumed  that  the  velocity  fluctuations  scale  with  the  velocity  jump, 
and  the  scale  with  the  FSL  width  (thickness);  and  it  was  also  assumed  that 
v  should  be  computed  on  the  dividing  streamline  (DSL).  Using  subscripts 
1  and  2  to  denote  the  "fast"  and  "slow"  sides  of  the  FSL  respectively, 
the  following  transition  Reynolds  number  based  on  FSL  thickness  (rather 
than  distance)  was  found: 


uihT  C'  Tdsl  k+1  A+1  T  % 

- TTMJ  '  T~i”  }  IT  -  f(Ml’  r0^ 


where  r^Mi)  expresses  the  dependence  of  the  fluctuations  on  Mi  and 


The  constant  C'  includes  the  magnitude  of  Re^o,  the  proportionality  constant 

between  A  and  the  FSL  thickness  h,  and  the  corresponding  constant  connecting 

u'  with  uj  -  u 2 •  Since  these  functional  connections  are  in  principle 

universal  for  each  self-preserving  flow,  the  dependence  of  h^  on  Mi,  A  and 

T  /T  can  be  then  found  from  eq.  (4),  while  its  precise  numerical  value 
02  0 1 

can  be  updated  as  more  information  on  the  ingredients  of  C'  becomes  avail¬ 
able  (C'  -  230  according  to  the  current  state  of  knowledge).  A  plot  of  eq. 
(4)  appears  on  Figure  2. 
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Reference  4  went  a  step  further  in  re-casting  eq.  (4)  in  a  way 
replacing  hy  by  the  "transition  distance"  Xy  from  the  origin  of  the 
laminar  FSL.  To  do  this,  it  was  necessary  to  find  a  relation  between 
h  and  X,  that  is  to  utilize  a  formula  for  the  FSL  growth  as  a  function 
of  distance  from  its  origin.  The  growth  formula  utilized  (in  Reference 
4)  is  plotted  on  Figure  3  and  was  due  to  Moeny  (Reference  3): 


4-  =  G(Mi’ 

x  o2  o, 


(6) 


where  Re  was  again  based  on  the  fast  side  properties.  Thus  eq.  (4)  gave 


Re 
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where  Z"  is  another  "universal"  constant  analogous  to  C'  and  simply  related 
to  it  (Reference  4).  If  C'  =  234,  C"  =  54,800.  Since  all  other  quantities 
on  the  r.h.s.  of  (7)  are  only  functions  of  Hu  \  and  T  /T  ,  the  transition 

02  Oi 

distance  Xy  from  the  origin  of  the  laminar  FSL  can  be  found  once  these 
three  parameters  of  the  flow  are  given.  Comparison  of  eq.  (")  with  data, 
made  in  Reference  3,  wa:.  satisfactory.  However,  there  are  important  pre¬ 
cautions  for  the  use  of  eq.  (7)  with  trail ing-edge  flows,  which  will  be 
explained  in  Section  2.4.  Eq,  (7)  is  plotted  on  Figure  4. 


Before  the  plans  for  further  verifying  eq.  (4)  or  eq.  (7),  which  led 
to  the  present  experiments,  are  described,  it  is  necessary  to  summarize  the 
reservations  which  a  fluid  dynamicist  should  have  in  using  the  present 
theory,  and  especially  equations  (4)  and  (7): 


a.  The  nature  of  these  formulas  as  estimates  of  minimum  transition 
Reynolds  numbers  (or  lengths)  should  be  kept  in  mind,  and  the  lack 
of  a  sufficient  condition  noted. 

b.  The  algebraic  forms  shown,  and  the  values  of  C'  and  C",  pre¬ 
suppose  fully  developed  turbulent  flows  following  transition. 

c.  The  data  available  for  establishing  relations  between  A  and  h. 
u'  and  ug  -  u(0)  etc.  are  extremely  meagre. 
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d.  The  values  of  C'  and  C"  are  flexible,  as  noted. 

e.  No  provision  is  made  for  lateral  variations  of  u',  A  and  the 
fluid  temperature  across  the  layer  height  or  width. 

f.  In  predicting  transition  events  or  comparing  the  equations  with 
data,  note  that  transition  to  turbulence  in  a  shear  flow  is  an 
event  protracted  over  several,  often  many,  shear  layer  thicknesses, 
usually  preceded  by  instabilities  and  followed  by  a  period  of 
turbulence  non-equilibrium.  In  the  present  method  this  zone  is 
contracted  to  a  point. 

Within  the  above  reservations  equations  (3),  (4)  and  (7)  represent 
transition  predictions  for  wakes  and  free  shear  layers.  In  the  next 
section  we  will  argue  that  the  mixing  flow  downstream  of  a  T.E.  is  actually 
a  combination  of  these  two  flows,  and  that  some  combination  of  the  equations 
will  be  needed  to  handle  such  a  geometry. 

2.4  Transition  in  a  FSL  Shed  From  a  Partition 


In  the  preceding  we  derived  formulas  for  FSL  transition  according  to 
which,  if  Mi  and  T02/T0l  are  fixed,  advance  knowledge  of  the  fast-side 
unit  Reynolds  number  and  the  speed  ratio  X  =  (ui-u2)/(ui+u2)  (or  r  =  u2/ui) 

will  determine  the  FSL  thickness  h  at  which  the  flow  will  become  turbulent. 
Still  within  the  framework  of  this  theory,  however,  it  is  clear  that  this 
situation  cannot  persist  when  X  approaches  zero,  if  the  FSL  is  shed  from 
a  partition.  Here  the  latter  term  signifies  the  non-equilibrium  shear 
layer  which  carries  along  the  wake  component.  Specifically,  the  theory 
(cf.  Figures  2,  3  etc.)  indicates  that  at  X  =  0  the  Re^y  becomes  infinite, 
as  appropriate  when  the  only  mechanism  for  shear  production  is  the  usual 
asymptotic  FSL  profile.  In  reality,  a  finite  RehT  limit  is  reached  at 
X  =  0  because  shear  is  produced  by  the  wake  component. 

This  limit  is  illustrated  on  Figure  5.  When  X  is  decreased  u2  approaches 
ui,  so  that  the  flow  profile  changes  fundamentally  to  one  of  a  wake.  Thus 
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at  \  =  0  the  profile  instability  becomes  wake  instability,  and  transition 
guidelines  must  be  drawn  from  eq.  (3)  rather  than  eq.  (4).  The  former 
equation  translates  into  a  straight  line  at  constant  Rehy  as  shown  on  the 
Figure,  whose  position  depends  on  the  wake  Reynolds  number  Re  and  which, 
in  turn,  depends  on  Re{  and  h;  the  dependence  of  Rehy  on  Rew  is  actually 
quite  complex  as  will  be  seen  shortly.  In  the  meantime  the  solid  line 
of  Figure  5,  joining  the  two  limiting  cases  becomes  the  complete  estimate 
for  transition  in  the  FSL  beyond  the  trailing  edge. 

To  calculate  specific  applications  of  Figure  5,  one  can  transform 
eq.  (3)  (especially  using  Reference  7  in  conjunction)  from  the  transition 
distance  Xy  into  a  transition  Reynolds  number  Re^y: 

RehT  =  (constant)  Mg2  ew _ 

0  •  8 

-  1) 

where 

=  .  Re,  •  U,C0h 

as  already  indicated.  The  "wake  drag"  is  defined  as 

Cph  =  2  x  (sum  of  T.E.  momentum  thicknesses) 

=  2(9:  +  02)  =  2  H  (10) 

The  "wake  Reynolds  number"  Re  is  thus  a  parameter  of  the  problem  while 

W 

ReQ  =  335  (Reference  7).  The  constant  appearing  in  eq.  (8)  has  the  value 
of  0.0912  if  in  eq.  (3)  we  put  y  =  1.4,  o  =  0.85,  Q  =  0,  and  k  =  0.75  for 
adiabatic  air  wakes.  With  these  values  Re^y/Me2  plots  as  in  Figure  6. 

We  shall  return  to  this  figure  momentarily;  first  we  shall  observe 
its  consequences  on  FSL  transition.  As  an  example,  Rehy  vs.  X  is  replotted 
from  Figure  2  on  Figure  7  for  the  case  M!  =  3,  using  C'  =  234.  At  X  =  1 , 
Re^y  is  about  7500  and  decreases  as  X  decreases.  When  X  is  very  small  the 
FSL  profile  progressively  approaches  the  wake  profile,  and  the  transition 
Reynolds  number  thus  approaches  a  finite  limit  given  by  eq.  (3)  or  its 


(8) 


(9) 


equivalent  eq.  (8)  (and  Figure  1).  This  limit  depends  on  Re  ,  and  its 

W 

lowest  possible  value  is  about  200QMer!  -  18,000  according  to  Figure  6. 

For  a  smaller  Rew,  say  400,  the  limit  increases  to  about  45,000.  These 
two  examples  are  drawn  in  Figure  7  as  straight  lines  parallel  to  the  X 
axis.  Thus  at  very  small  X  the  FSL  transition  behavior  commences  at 
this  "wake-controlled  limit”,  and  as  X  increases  toward  unity  the  behavior 
changes  to  conform  to  the  "FSL-contrcl led"  curve  originating  in  eq.  (3), 
Clearly,  these  two  extremes  must  be  joined  smoothly  by  a  path  such  as 
shown  by  dashed  lines  on  Figure  7. 


Further  insight  (and  agreement  with  intuition)  comes  from  further 
study  of  Figures  6  and  7.  To  begin  with,  since  Rew  =  4Re5j£  (where  Re^ 
is  the  momentum  Reynolds  number  of  the  boundary  layer  at  the  T.E.)  one 
notes  that  Re  >  2000  or  so  on  Figure  6  corresponds  to  Re,,T.-  >  500.  In 
most  laboratory  and  practical  applications  this  Rea  level  marks  the  onset 
of  turbulence  in  the  boundary  layer  itself.  Thus  Re  >  2000,  approximate 
has  no  corresponding  laminar  wake  flow,  and  that  portion  of  the  curve 
on  Figure  6  is  of  no  physical  significance.  In  the  portion  Re  <  2000 

W  v 

one  can  divide  the  flow  roughly  into  three  categories: 


-For  Re^  <  300  no  turbulent  wake  flow  is  possible  (Reference  7). 

Thus  no  limit  to  the  FSL-controlled  transition  behavior  will  appear 
on  Figure  7,  at  any  Me.  Here  the  wake  component  is  so  feeble 
that  it  has  no  opportunity  to  destabilize  and  trip  the  flow. 


--For  intermediate  Re,.,  say  400  <  Rew  <  800,  a  typical  limit  is 

w  w 

shown  on  Figure  7.  The  crossing  occurs  at  X  =  0.02. 

--For  large  allowed  Rew,  we  say  1000  <  Rew  <  2000,  the  RehT  attains  a 
minimum  depending  on  Mach  number.  For  Me  =  3  Figure  7  shows  this 
minimum  (at  18,000)  to  be  still  quite  larger  than  the  X  =  1  limit 
of  the  FSL-controlled  transition  curve.  The  crossing  occurs  at 
X  =  0.07. 


The  present  purpose  is  to  create  conditions  in  the  laboratory  by 
which  the  predictions  detailed  above,  especially  those  in  this  section, 
can  be  tested.  The  planning  of  experiments  for  tnis  purpose  is  explained 
in  the  next  section. 


3 .  The  Experiments:  Motivation  and  Design 

The  purpose  of  the  experimental  work  described  here  was  to  detect 
and  correlate  laminar-turbulent  transition  events  in  the  FSL  separating 
two  parallel  streams  issuing  from  adjacent  nozzles.  In  order  to  reduce 
the  scope  of  the  program  to  a  manageable  size  the  flows  chosen  were 
steady,  isothermal  and  homogeneous;  to  provide  data  relevant  to  fluid 
laser  design  the  flows  were  supersonic;  and  to  enable  comparison  with 
the  transition  predictions  noted  above  (and  the  laminar  FSL  theory  shown 
in  the  Appendix)  the  trailing  edge  of  the  partition  separating  the 
streams  was  made  sharp. 

Despite  its  shortcomings  listed  in  the  previous  Section,  sufficient 
resemblance  between  the  few  previous  available  data  and  the  "transition" 
theory  warranted  the  use  of  the  latter  as  a  design  tool  for  the  experi¬ 
ments.  Figures  8  and  9,  reproduced  from  the  original  proposal  for  this 
work,  can  be  used  to  illustrate  the  following  discussion.  Figure  8  is 
a  qualitative  picture  of  the  events  expected  to  occur  when  two  streams 
of  Mach  numbers  Mi  and  M2  (M2  <  Mi)  are  made  to  flow  parallel  to  each 
other  after  issuing  from  the  adjacent  nozzles.  To  simplify  the  picture 
it  is  first  assumed  that  the  transition  movement  in  the  boundary  layer 
of  the  top  surface  of  the  partition  coincides  with  the  corresponding 
movement  of  the  layer  on  the  bottom  surface;  thus  there  exists  a  unit 
Reynolds  number  Ref  below  which  the  entire  flow  coming  off  the  T.E.  is 
laminar.  Changes  in  Re'  should  then  affect  transition  as  follows:  at 
very  large  Re'  the  nozzle  wall  boundary  layers  are  turbulent,  and  as 
Re'  is  lowered  transition  moves  downstream  along  the  nozzle  surfaces 
till  it  reaches  the  T.E.  The  particular  shape  of  this  variation  (curve 
"A"  in  Figure  8)  is  not  a  problem  addressed  here,  although  it  should 
be  of  interest  for  applications.  Curve  "A"  can  be  assumed  known  for  a 
specific  nozzle  geometry  and  flow  conditions. 

When  on  the  other  extreme  Re'  is  very  low  eq.(7)  implies  that  the 
FSL  remains  laminar  at  far  distances  downstream  of  the  T.E.  In  this  case 
the  FSL  reaches  the  equilibrium  laminar  profile  discussed  in  Section  2 
and  the  Appendix;  the  "memory"  of  the  wake  component  is  lost  both  in  the 
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pre-transitional  laminar  and  the  post-transitional  turbulent  flow.  In 
this  case  eq.  (7)  is  supposed  to  apply  and  its  Ref  -  XT  prediction  is 
shown  as  curve  "C"  in  Figure  8. 

The  problem  lies  in  connecting  curves  A  and  C.  If  eq.  (7)  was  valid 
near  the  T.E.  (ruled  as  unlikely  by  the  preceding  comments)  a  decrease 
in  Re{  would  cause  transition  to  "jump"  from  the  T.E.  to  curve  C  and  then 
follow  the  latter  as  Re{  is  further  decreased.  The  invalidity  of  eq.  (7) 
near  the  T.E.  implies,  instead,  that  a  more  gradual  joining  of  A  and  C 
would  occur,  hypothesized  by  curve  "B"  on  Figure  8.  The  precise  track 
of  this  curve  was  the  objective  of  this  work. 

The  path  of  curve  B  was  of  course  unknown  other  than  that  it  should 
start  at  the  T.E.  and  asymptotically  approach  curve  C  if  eq.  (7)  was 
valid.  The  design  of  the  experiment  required  fixing  M2  so  that  the  curve 
B  was  most  easily  detectable.  Adopting  for  the  moment  the  predicted 
transition  distance  X-j.  of  eq.  (7)  as  valid,  computations  were  made  of 
the  transition  occurence  in  the  MSU  Supersonic  Wind-Tunnel  (SWT)  taking 
one  stream  to  be  at  Mi  =  3  and  the  other  at  a  series  of  lower  Mach  numbers 
M2  <  Mi.  The  results  consist  of  curves  in  the  Ref  -  X  plane,  each  of 
which  separates  the  turbulent  flow  above  it  from  laminar  flow  below  it, 
as  shown  on  Figure  9.  The  same  figure  shows  observability  "window"  (or 
"test  aperture")  for  the  FSL,  i.e.  it  sets  limits  on  the  upper  and  lower 
possible  Ref  (because  of  PQ  limitations)  and  on  the  maximum  downstream 
distance  X  because  of  the  tunnel  size. 

A  serious  design  difficulty  arose  in  establishing  the  virtual  origin 
of  ie  laminar  FSL,  on  which  XT  of  eq.  (7)  is  based.  Since  no  character¬ 
istic  scale  length  appears  in  the  formulation  of  the  equilibrium,  self¬ 
similar  FSL  this  issue  was  not  important  in  earlier  comparisons  of  the 
theory  with  data  (cf.  Reference  3,  Figure  1)  because  the  latter  involved 
mainly  FSL's  generated  by  crossing  shock  waves  (References  8  and  9).  In 
planning  tests  with  actual  FSL's  beginning  with  finite  thickness  at  some 
T.E.,  one  can  at  most  assume  that  the  flow  began  at  some  fictitious 
(virtual)  origin  upstream  of  the  T.E.  As  a  last  resort,  this  unsatis¬ 
factory  process  was  used  in  the  design  (Reference  10).  The  estimated 
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virtual  origin  location  is  shown  on  Figure  9,  marking  the  location 
Xj  =  0.  The  method  of  computing  the  virtual  origin,  explained  in 
Reference  10,  will  be  contested  by  the  data  presented  later  in  this 
report.  In  the  meantime,  since  meaningful  measurements  could  be 
made  only  for  M2  curves  near  the  center  of  the  test  aperture.  Figure 
9  indicated  the  choice  of  M:  =  1.5. 


4 .  Matrix  of  Experimental  Co nd vt jo ns 

As  a  result  of  the  transition  expectations  of  Section  2  and  the 
design  considerations  expressed  in  Section  3,  the  following  initial 
experiment  conditions  were  set: 

M,  =  3 
M;  =  1.5 

p0/p0, =  10 
To  =  100°F  =  V  To2 

=  300  -  600  mm  Hg.  abs. 

At  these  conditions  the  expected  X  was  0.18,  the  expected  ReXT/C"  =  24 
(ReXT  =  1,300,000)  according  to  Figure  4  and  Re^/C'  =  45  (RehT  =  10,500) 
according  to  Figure  2.  Furthermore,  Figure  9  indicated,  at  \  =  0.18. 
that  transition  would  be  controlled  by  the  FSL  character  (profile)  of  the 
flow  rather  than  by  its  wake  character.  Finally,  according  to  the  expec¬ 
tations  of  Figure  7,  by  changing  pQ^  from  600  to  300  mm.  Hg. ,  the  transition 
point  should  be  observed  to  move  from  the  T.E.  to  about  12"  (30  cm.)  down¬ 
stream. 

The  fixed  ratio  P0i/P0  =  10  indicated  that  the  static  pressures  at 
the  exit  of  the  two  nozzles  would  be  equal,  so  that  the  FSL  was  expected 
to  progress  down  the  middle  of  the  tunnel  undeflected. 

We  shall  see  further  below  that  a  variety  of  reasons  combined  to  change 
the  matrix,  especially  the  levels  of  M2 ,  p  /p  and  X. 
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5.1  Wind-Tunnel 


These  experiments  were  conducted  in  the  Supersonic  Wind-Tunnel  at 
Montana  State  University  (MSU/SWT).  The  facility  is  described  in  detail 
in  Reference  11.  A  summary  of  its  salient  features  is  given  below. 

The  wind-tunnel  has  been  in  almost  uninterrupted  use  since  the  1960's. 

It  consists  of  a  continuous,  open-circuit  flow  facility  with  the  components 
upstream  of  the  test  section  consisting  of  the  atmospheric  inlet,  a  silica- 
gel  air  dryer  with  a  dewpoint  capability  of  -35'F  doubling  as  an  air  heater, 
a  throttling  valve  controlling  the  supply  (stagnation)  pressure  and  a  large 
stilling  tank.  Downstream  of  the  test  section  the  air  passes  through  super¬ 
sonic  and  subsonic  diffusers  and  the  pumping  plant  driving  the  tunnel  and 
is  thence  discharged  to  the  atmosphere.  A  schematic  of  the  tunnel  circuit 
is  shown  on  Figure  10. 

In  its  normal  configuration  the  test  section  is  composed  of  a  two- 
dimensional  DeLaval  nozzle  designed  for  exit  M  =  3.0,  terminating  in  a 
3.1"  x  3.2"  (7.87  x  8.13  cm.)  rectangular  test  section.  The  removable 
sidewalls  are  made  of  single  pieces  of  optical -qual i ty  glass  stretching 
from  upstream  of  the  nozzle  throat  to  the  diffuser  entrance.  Diagnostic 
probes  can  be  suspended  on  vertical  struts  passing  through  the  test-section 
ceiling  and  slaved  to  the  electromechanical  actuator  system  (Section  5.3.6). 
The  test  section  flow  is  very  uniform  within  r0.04  of  M  =  3. 

Control  of  the  tunnel  is  maintained  from  a  console  with  pQ  set  point 
to  ±1  mm  Hg.  abs.,  Tq  ret  point  to  ±1°F,  and  dewpoint  readouts  for  points 
across  the  air  dryer.  Maximum  pQ  is  atmospheric  (usually  625  mm  Hg.)  and 
maximum  Tq  of  order  150°F.  Maximum  and  minimum  unit  Reynolds  attained 
depend  on  configuration  but  vary  from  80,900  cm"1  to  20,000  cm"1  respectively. 

5 . 2  Modification  of  the  SWT  for  the  Present  Measurements 

In  general  outline  the  modification  made  to  the  SWT  test  section  for 
the  present  measurements  is  shown  on  Figure  11,  and  is  also  illustrated  in 
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the  photographs  of  Figures  12  and  13.  Basically  the  Mach  3  nozzle  is 
"split  in  half"  by  positioning  a  flat  splitter  plate  with  a  sharp  T.E. 
along  its  centerplane.  A  "lower"  channel  is  thus  formed,  visible  in 
Figures  11,  12  and  13  of  half  of  the  original  Mach  3  nozzle,  with  its 
floor  being  the  lower  block  of  the  latter,  and  the  flat  plate  forming  its 
ceiling.  The  characteristics  method  predicts  normal  flow  for  this  type 
of  configuration,  as  already  verified  in  the  past  for  this  tunnel 
(Reference  12).  This  "lower"  channel  is  here  frequently  called  the  "fast" 
or  "Mach  3"  channel  herein  and  its  flow  properties  are  subscribed  by  "1". 

The  upper  surface  of  the  flat  plate  was  contoured  to  produce  a  Mach 
number  M2  <  Mi  in  the  "upper"  channel  shown  on  Figures  11,  12  etc.,  whose 
ceiling  was  formed  by  another  flat  plate  parallel  to  the  SWT  principal 
axis;  this  is  hereafter  called  the  "slow"  channel  and  its  properties  are 
subscribed  by  "2".  Air  to  provide  this  channel  with  flow  was  diverted 
from  the  main  stilling  tank,  with  a  manually-controlled  valve  interposed, 
as  shown  on  Figure  11.  This  valve  controlled  pQ  ,  the  low-speed  channel 
stagnation  pressure,  while  the  main  SWT  throttling  valve  mentioned  previ¬ 
ously  controlled  po^  and  pQ  together.  By  means  of  this  arrangement  pQi 
and  pQ  could  be  set  independently;  or  pQ^  and  pQ  could  be  raised  or 
lowered  in  unison  while  their  ratio  remained  fixed  at  a  desired  value. 

The  stagnation  temperature  Tq  was  always  the  same  for  both  channels.  In 
all  cases  the  splitter  plate  spanned  the  test  section  completely. 

The  important  parameters  Mi,  M2,  9i,  92,  Ref ,  Re2  of  the  measurements 
were  set  by  the  upper-side  contour  of  the  splitter  plate  and  the  distance 
xi  between  its  T.E.  and  the  Mi  =  3  throat.  The  various  configurations 
(or  "modes")  resulting  from  combinations  of  these  are  shown  on  Table  I. 

The  progression  from  mode  I  to  mode  III  reflects  the  passage  of  calendar 
time  and  the  desire  to  improve  the  M2  nozzle  contour  and  to  solve  the  SWT 
choking  problems;  details  on  the  flows  obtained  with  each  mode  will  be 
given  below.  Contour  coordinates  for  the  M2  nozzle  were  provided  by  the 
short-nozzle  theory  of  Shames  and  Seashore  (Reference  13).  The  nozzle 
coordinates  for  modes  II  and  III  are  shown  on  Tables  II  and  III.  The 
trailing  edges  were  made  as  sharp  as  possible.  Figures  11  and  12  show 
the  geometries  of  the  modes  listed  on  Table  I. 
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TABLE  II.  Coordinates  of 
the  M=1 .6  Nozzle  (Modes  I,  II) 


TABLE  III-  Coordinates 

of  the  M=2.25  Nozzle  (Mode  III) 


x  (inch) 


0 

.6805 
.8880 
1.0545 
1.1526 
1.2547 
1 .4253 
1 . 5865 
1.6775 
1.8853 
2.0942 
2.474 
2.8414 
3.2729 
3.7117 


y  (inch) 


.35 

.2571 

.2298 

.2081 

.1954 

.1822 

.1607 

.1411 

.1303 

.1069 

.0852 

.0512 

.0263 

.0076 

0 


x  (inch) 


0 

.52219 
.640224 
.709047 
.757126 
.795076 
.82692 
.854658 
.901887 
.94183 
.976921 
1.02331 
1 .06443 
1.1018 
1.13635 
1.1687 
1.19929 
1.28917 
1.34572 
1.45038 
1.49977 
1.68586 
1.86225 
2.05676 
2.27339 
2.49424 
2.7219 
3.08047 
3.59892 
4.1743 
4.81843 


y  (inch) 


.87864 

.72387 

.68893 

.668604 

.65443 

.64327 

.63392 

.62579 

.61199 

.60036 

.59018 

.5768 

.565 

.554345 

.54456 

.53544 

.52688 

.50203 

.48666 

.4588 

.4459 

.39912 

.35739 

.31441 

.2703 

.22934 

.1912 

.13919 

.08014 

.0349 

.0068 


During  the  tests  of  the  modes  I,  II,  II  ,  and  III,  various  diffuser 

d 

combinations  were  attempted,  including  a  variable  one  consisting  of  two 
plates  hinged  at  their  junction,  so  that  a  variable  wedge  geometry  could 
be  generated  by  pulling  wires  through  the  test  section  ceiling. 

5.3  Instrumentation 

5.3.1  Schlieren  Optical  System 

The  SWT  Schlieren  system  was  extensively  used  to  discern  the  overall 
flow  features  including  choking  phenomena,  detect  transition  and  even 
confirm  the  existence  of  vortical  instabilities. 

The  system  is  based  on  the  light  from  a  tungsten  source  passing  in 
parallel  fashion  through  the  test  section  windows.  Photography  of  the 
image  was  done  both  in  continuous  exposure  and  by  spark;  the  latter  has 
an  equivalent  exposure  time  of  one  microsecond.  Schlieren  observations, 
which  will  be  later  shown  by  photographs,  were  greatly  aided  by  the  wide 
field  of  view  of  the  flow  made  possible  by  the  large  window  area  of  the 
test  section. 

5.3.2  Hot-film  Anemometer 

A  hot-film  anemometer  was  used  to  measure  the  location  of  transition 
and  to  make  qualitative  measurements  of  the  turbulent  intensity  and  its 
spectrum.  The  film  consisted  of  platinum,  approximately  0.002"  x  0.02" 
(0.005  x  0.05  cm)  in  dimensions,  and  was  deposited  on  the  leading  edge  of 
a  double  wedge  ground  on  the  tip  of  a  0.02"  (0.05  cm)  diameter  kimax  glass 
rod.  The  probe  was  connected  to  either  a  Transmetrics  6401  or  an  ADP1213 
hot-wire  heating  circuit  and  amplifier  combinations.  These  circuits  act 
both  as  generators  of  the  heating  current  (usually  about  20  ma  d.c.)  of 
the  probe,  and  also  as  amplifiers  of  its  a.c.  signal.  The  upper  frequency 
response  of  these  amplifiers  were  320  and  3,000  KHz  respectively,  with 
reponse  compensation  as  high  as  100  msec  and  gain  to  50.000.  For  many  of 
the  observations  the  film  signals  were  first  recorded  at  120  ips.  on  a  7- 
channel  tape  recorder  with  FM  wideband  response  to  400  KHz. 
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Only  a  fraction  of  this  capability  was  utilized,  since  no  quantitative 
turbulence  measurements  were  intended.  Transition  data  were  taken  by 
judging  the  frequency  extent  of  the  signals  and  by  their  intensity  within 
their  active  frequency  band.  Oscilloscopes,  r.m.s.  volt-meters  and  wave 
analyzers  were  used  to  make  this  judgement.  Especially  useful  in  this 
respect  was  the  Rockland  FFT512/S  fast  Fourier  transform  computer  which 
provided  smoothed  spectra  by  repeated  averaging. 

5.3.3  Static  Pressure  Probe 

Static  pressure  measurements  were  made  with  a  probe  consisting  of  a 
0.056  cm-dia.  tube  with  3  holes  drilled  around  its  circumference  at  a  distance 
of  0.46  cm  from  its  sharpened  tip.  Signal-processing  gear  for  this  probe 
were  identical  to  that  connected  to  the  pitot  probe,  and  will  be  described 
below.  The  minimum  lateral  resolution  of  this  probe  equalled  its  diameter. 

No  static  pressure  measurements  were  taken  in  the  region  downstream  of  the 
T.E.  and  for  a  length  of  0.5  cm.  from  it. 

5.3.4  Pitot  Probe 

The  pitot  probe  was  the  principal  instrument  for  the  measurement  of 
mean  (average)  flow  quantities.  It  consisted  of  a  0.008"  (0.02  cm)  diameter 
tube,  flattened  and  chemically  etched  at  the  tip  to  a  lateral  (y)  dimension 
of  0.005"  (0.013  cm.).  Since  the  shear  layer  width  was  of  order  0.15"  (0.38 
cm.)  in  width,  this  afforded  a  lateral  resolution  of  0.15  0.005  =  30 

which  is  considered  quite  adequate.  The  probe  tube  was  telescoped  into 
tubes  of  progressively  larger  diameter,  for  rigidity,  and  terminated  in  a 
bullet-shaped  housing  which  encapsulated  a  Kulite  XTH-1-190-5A  pressure 
transducer.  The  transducer  housing  itself  was  immersed  in  the  flow,  thus 
shortening  the  distance  between  probe  tip  and  transducer  and  thereby  greatly 
decreasing  the  response  time. 

The  transducer  was  energized  by  15.7  volts  d.c.  and  had  a  pre-ampl ified 
sensitivity  of  22.4  mv/v/mm.  During  data  recording  the  transducer  output 
was  amplified  by  100  and  dynamically  damped  to  eliminate  electronic  noise 


and  drift.  The  transducer  was  calibrated  at  least  once  every  day  during 
the  tests. 


5.3.5  Total  Temperature  Probe 

The  total  temperature  probe  consisted  of  a  type  K  thermocouple  bead, 
0.013  cm.  in  diameter,  connected  to  a  digital  temperature  indicator  outside 
the  flow.  The  calibration  of  this  instrument  v.  Reynolds  number  will  be 
given  during  the  discussion  of  the  data. 

5.3.6  Electromechanical  Positioning  and  Actuating 

All  flow- immersible  probes  were  held  with  their  tip  on  the  vertical 
mid-span  plane  of  the  flow.  They  were  supported  in  two  degrees  of  freedom 
by  the  SWT  electromechanical  actuator,  consisting  of  a  horizontal  hollow 
strut  held  in  place  by  two  vertical  struts  protruding  through  the  ceiling. 
The  vertical  struts  are  actuated  in  the  vertical  (y)  direction,  while  an 
independent  motor  drove  the  hollow  strut  along  the  flow  (the  x  direction). 
The  actuator  control  provided  start/stop,  speed  and  direction  commands  for 
either  of  these  two  motions,  as  well  as  digital  counter  and  d.c.  voltage 
position  indications. 

For  the  mean  flow  measurements  the  data  were  obtained  in  the  form  of 
"profiles",  i.e.  traverses  along  y  at  constant  x.  The  SWT  electromechanical 
actuator  system  allowed  this  traverse  to  be  done  fully  automatically.  The 
vertical  motion  is  geared  to  a  flywheel  with  adjustable  spokes  protruding 
from  its  periphery  which  periodically  interrupt  a  photodetector  beam,  each 
interruption  transmitting  an  electric  "read"  command  to  the  SWT  Data  System. 
By  adjusting  the  angular  separation  of  the  spokes  around  the  flywheel,  a 
"read"  command  at  each  of  the  vertical  steps  mentioned  above  is  obtained. 

In  addition,  and  to  ensure  adequate  response  time  for  the  sensor,  the  system 
automatically  produces  a  pause  between  sensor  arrival  at  a  new  vertical 
point  and  the  "read"  command  at  that  point,  the  pause  being  typically  15 
seconds  in  length.  The  probe  signals  were  acquired  and  stored  by  the  SWT 
Data  System  which  basically  consists  of  a  Spectral  Dynamics  Model  SD-133 
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8-channel  A/D  converter  and  a  Texas  Instruments  Silent  790  ASR  computer 
terminal  with  dual  cassette  storage.  This  combination  allows  storage  on 
cassette  as  well  as  real-time  hard-copy  printing  of  the  line  prompted  by 
each  "read"  command.  Alternatively  an  Intertek  Superbrain  microcomputer 
with  dual  5-inch  disk  storage  was  also  used  in  series,  to  provide  CRT 
display  of  the  data  and  disk  storage. 
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TABLE  I 


NOZZLE  CONFIGURATION  MODES 


I 

II 

I  la 

III 

Mi  nominal 

3 

3 

2.7 

2.88 

Mi  actual 

3 

3 

- 

2.77 

Type  1  nozzle 

contoured3 

contoured3 

contoured3 

contoured3 

M2  nominal 

1.56 

1.6 

1.6 

2.25 

M2  actual 

1.47 

- 

1.74 

2.18 

Type  2  nozzle 

straight 

contoured1 

contoured2 

contoured 3 

Xi,  inch  (cm) 

10.4"  (26.4) 

10"  (25.4) 

5"  (12.7) 

6.65  (16.9) 

x2,  inch  (cm) 

4.15"  (10.5) 

3.71"  (9.42) 

3.71"  (9.42) 

4.8  (12.2) 

X  nominal 

0.167 

0.159 

0.139 

.0537 

X  actual 

0.187 

.0542 

1  Template 

2  Template  modified  by  hand 

3  Numerically-controlled  mill . 
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6.  Results 


6. 1  Preliminary  Results  With  Modes  I  and  II 

The  initial  measurements,  as  per  the  planned  matrix  of  Table  I, 
were  carried  out  with  Modes  I  and  II  for  which  the  lower-speed  channel 
used  a  Mach  1.6  nozzle.  This  was  primarily  a  "shakedown"  experiment,  aimed 
at  judging  the  overall  performance  of  the  modified  SWT  test  section.  As 
such,  much  effort  went  into  observing  design  defects,  into  making  modifica¬ 
tions  and  into  explaining  the  overall  flow  phenomena  observed. 

Modes  I  and  II  aimed  at  achieving  the  flow  conditions  shown  on  Figure  9, 
i.e.,  a  long  FSL  parallel  to  the  tunnel  axis,  with  Mi  =  3,  Ma  =  1.6  and 
PQ  /PQ,  =  10’  1°  an  overall  sense  all  the  hardware  performed  properly; 

observations  made  with  Mode  I  (see  Table  I)  as  well  as  with  Mode  II,  shown 
in  structural  outline  on  Figure  13,  produced  the  type  of  flow  seen  on 
Figure  14.  However  a  number  of  difficulties  were  encountered  at  the  outset. 

First,  the  desired  P  /P  =10  condition  could  not  be  achieved  without 

oi  o  2 

flow  breakdown  at  the  nozzle  exits,  because  of  inadequate  diffuser 
performance.  Second,  while  satisfactory  flow  was  achieved  at  off-design 
conditions  of  P  /P  <  10,  this  meant  an  increase  in  P  which  in  turn 

O 1  02  02 

caused  a  downward  deflection  of  the  splitter  plate;  in  turn,  this  caused 
an  increase  in  M;  and  a  decrease  in  Mi.  Third,  the  contour  of  the  M2 
nozzle,  fabricated  from  a  template,  proved  unsatisfactory  for  uniform  exit 
flow.  The  combined  result  of  these  problems  are  vividly  seen  on  Figure  14, 
in  the  form  of  a  deflected  FSL  which  in  addition  terminates  into  a  maze  of 
spurious  shock  waves  and  separated  wall  boundary  layers  before  the  diffuser 
entrance. 

Original  estimates  had  not  indicated  any  potential  problems  with  the 
SWT  diffuser.  This  had  been  encouraging  since  it  indicated  no  need  to 
modify  the  entire  diffuser  section  according  to  most  recent  "long-diffuser" 
ideas  for  dissimi lar-stream  flows  (e.g.  Reference  14).  Prolonged  efforts 
to  improve  the  diffusion  for  Modes  I  and  II,  including  the  hinged  variable 
diffuser  already  mentioned  in  Section  5.2,  produced  distinct,  but  not 
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drastic,  improvements.  It  should  be  kept  in  mind,  too,  that  the 
immersion  of  diagnostic  probes  in  the  test  section  (some  are  visible  in 
Figures  12  and  13)  provide  additional  obstacles  to  proper  diffusion. 

Relief  to  the  choking  difficulty  appeared  from  an  unexpected  source, 
which  will  be  described  briefly  here.  It  was  found  by  measurements  with 
the  pilot  tube,  that  the  trailing-edge  (initial)  FSL  thickness  was 
considerably  higher  than  originally  estimated  from  Reference  4  and  included 
in  Reference  10.  Specifically,  Figure  15  shows  the  T.E.  shear  thickness 
6j  +  5,  as  calculated  in  the  experiment  design  phase,  and  as  measured  for 
the  Mode  I - 1 1  tests.  These  measurements  show  considerable  scatter 
because  they  were  preliminary,  but  in  anyway  indicate  a  T.E.  FSL  thickness 
2-3  times  greater  than  anticipated.  Since  this  finding  had  important  effects 
on  the  experimental  plan,  an  attempt  to  support  it  by  independent 
theoretical  means  was  sought.  It  had  been  previously  found  that,  despite 
the  curvature  and  pressure  gradients  present  in  DeLaval  nozzles,  reasonable 
estimates  (to  less  than  15%)  of  their  boundary-layer  thickness  could  be 
made  using  standard  laminar-flow  flat-plate  theories  such  as  that  due  to 
Low  (Reference  15).  Figure  16  shows  that  Low's  theory  supported  data 
obtained  in  typical  nozzles  (in  this  case  the  MSU/SWT  regular  M  =  3  nozzle) 
and  verifies  that  the  original  FSL  estimates  were  inappropriate. 

The  reason  why  these  findings  were  important  to  the  diffuser 
performance  was  that  they  re-directed  the  search  for  transition  closer  to 
the  T.E.  (since  an  originally  thicker  FSL  generally  requires  a  smaller 
wetted  length  to  transition).  This  made  the  requirement  for  a  long  FSL 
unnecessary.  Flow  breakdown  phenomena  near  the  diffuser  entrance  could  be 
tolerated,  therefore,  so  long  as  a  few  inches  (or  even  a  few  centimeters) 
of  FSL  length  were  obtained  downstream  of  the  T.E.  Further  work  on 
improving  diffusion  was  therefore  stopped,  and  for  the  remainder  of  this 
report  we  deal  with  FSL's  of  the  type  and  extent  shown  on  Figure  14. 

The  setback  caused  by  the  findings  of  Figure  15,  however,  was  that 
they  negated  eqs.  (6)  and  (7);  thus  although  transition  predictions  on  the 
basis  of  thickness  (eq.  (4)  and  Figures  2  and  7,  for  example)  remained 


23 


valid,  predictions  based  on  wetted  length  (eq.  (7)  and  Figure  4)  did  not. 

To  continue  expressing  transition  expectations  on  a  Re'  -  x  plane,  one 
needed  a  valid  h(x)  relation  to  replace  eq.  (6).  According  to  the 
Appendix,  such  a  relation  can  be  obtained  from  the  exact  theory,  but 
requires  tedious  computations  for  each  given  parameter  group  (Mi,  etc.). 

Therefore,  it  was  decided  to  reconsider  the  tests  underway  in  terms  of 
the  Poi,  h  plots  shown  in  this  report  beginning  with  Figure  17.  In  this 
representation  it  is  recognized  that  Poi  and  Po?  may  be  arranged  so  that 
the  FSL  at  the  T.E.  will  deflect  up  or  down  by  an  angle  a.  Thus  any  set 
of  nozzle-exit  M  ,  M  ,  P  ,  P  etc.  produces  a  different  set  of  M, ,  M. 
etc.  where  the  latter  is  now  the  true  properties  on  either  side  of  the 
deflected  FSL,  as  is  shown  on  Figure  18.  In  turn,  this  allows  the 
computation  of  X  and,  from  Figure  2,  the  transition  Reynolds  number  RehT 
and  also  the  transition  thickness  h  since  Re‘  is  now  known.  In  Figure  17 
the  transition  thickness  h  is  thus  plotted  as  a  function  of  P  and  P 

?  i  o ; 

for  the  Mode  I- 1 1  nozzles,  and  the  same  is  done  on  Figure  19  for  the 

Mode  I  la  nozzle.  One  can  of  course  also  plot  on  such  graphs  the  corresponding 

transition  thickness  if  the  flow  "acted  like  a  wake"  (Section  2.4),  as 

well  as  the  initial  (trailing  edge)  h  =  Si  +  52  for  orientation  purposes. 

6.2  Instability  and  Transition  Results  for  Modes  I  and  II 

The  conclusions  from  the  work  reported  in  the  last  section  were 
that  the  observed  FSL's  of  several  inches  in  length,  were  suitable  for 
making  transition  observations;  the  mechanical  difficulties  mentioned, 
however,  limited  any  further  experimentation  with  this  mode  valuable  only 
in  an  exploratory  sense. 

The  major  flow  feature  explored  was  the  turbulence  content  of  this 
FSL  at  P  =  570  mm  Hg,  P  =  296  mm  Hg.  At  these  conditions,  chosen 
because  of  the  flow  steadiness,  the  flow  deflected  toward  the  fast-flow 
channel  (positive  a  in  Figure  18),  thereby  decreasing  the  speed  difference 
between  the  two  channels;  it  is  estimated  that  Mi  =  2.8,  M2  =  1.8  and 
thus  X  =  0.11.  It  is  also  estimated  that  both  boundary  layers  were  laminar 
at  the  T.E. 
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Figure  20  shows  profiles  of  the  hot-film  wideband  a.c.  output  across 
the  shear  layer,  for  progressively  increasing  distances  x  from  T.E. 

These  exhibit  two  distinct  “peaks"  on  the  sides  of  the  FSL ,  which  grow 
quickly  from  the  T.E.  to  a  distance  of  about  2.5  cm.  from  the  T.E.  and 
thereafter  decay.  Spectra  of  the  probe  signals  along  both  the  top  and 
the  bottom  "peaks"  are  shown  on  Figures  21  and  22.  The  first  trace  on 
each  figure  is  taken  so  close  to  the  T.E.  that  they  can  be  considered  to 
be  at  x  =  0,  i.e.  in  the  two  T.E.  boundary  layers.  These  spectra  shed 
additional  light  into  what  causes  the  wideband  signals.  It  is  seen  from 
Figure  21  that  an  "energy  peak"  appears  within  a  fraction  of  an  inch  past 
the  T.E.,  which  centers  at  about  75  KHZ.  This  peak  actually  appears  first 
on  the  upper  (low-speed)  side  of  the  FSL,  and  is  accompanied  also  by 
signal  increases  at  the  higher  frequencies.  Note  that  the  frequency  of 
the  peak  does  not  change  as  far  from  the  T.E.  as  3.6"  (9.1  cm).  By  that 
time  it  is  seen  from  the  spectra,  that  the  oscillation  has  quite  well 
vanished  from  the  upper  (low-speed)  side  of  the  FSL,  while  it  persists  in 
the  lower  (high-speed)  side. 

The  following  conclusions  were  drawn  from  this  measurement: 

(a)  The  transition  picture  is  complex  and  involves  asymmetries  and 
non-uniformities  both  in  the  spectral  and  spatial  sense. 

(b)  If  turbulence  is  identified  as  a  spectral  activity  at  high 
frequencies,  it  appears  that  transition  sets  in  quite  abruptly 
within  one  centimeter  or  so  from  the  T.E. 

(c)  An  intense  periodic  phenomenon  appears  as  soon  as  the  two  flows 
meet  just  beyond  the  T.E.  The  periodicity  decays  quickly  on  the 
slow  side  of  the  FSL  but  remains  practically  intact  for  long 
distances  on  the  fast  side  of  the  FSL. 

These  conclusions  led  to  a  tentative  picture  of  events  in  the  FSL, 
which  consists  of  the  FSL  forming  initially  as  a  laminar  flow,  but  almost 
immediately  changing  to  one  which  is  both  turbulent  and  oscillatory.  A 
very  close  parallel  in  low-speed  flows  can  be  found  in  Reference  16. 


As  for  transition.  Figure  17  had  already  predicted  that  for  this 
mode  the  FSL  forming  at  the  T.E.  was  already  of  a  thickness  which  would 
require  it  to  be  turbulent,  according  to  the  asymptotic  theory.  The  flow 
should  thus  become  turbulent  as  soon  as  possible  beyond  the  T.E.,  which  is 
confirmed  by  the  above  observations. 

6.3  Results  With  Mode  III 

6.3.1  General 


This  configuration,  shown  on  Figure  13  and  Table  III,  was  designed  to 
produce  initial  boundary  layers  at  the  T.E.  thinner  than  would  be  required 
for  transition.  With  the  latter  thus  expected  some  distance  downstream  of 
the  T.E.  it  was  thought  that  a  clearer  picture  would  be  obtained  of  the  h^ 
(Re')  dependence,  for  which  predictions  appear  on  Figure  23. 

Figure  24  shows  Schlieren  photos  of  typical  flows  obtained  with  this 
configuration.  Absence  of  shock  waves  and  other  disturbances  indicates 
proper  nozzle  flow. 

Study  of  Mode  III  began  with  an  effort  to  see  if  the  boundary  layers 
merging  at  the  T.E.  were  laminar  or  turbulent;  since  the  present  program 
requires  initially  laminar  FSL  at  the  T.E.,  it  follows  that  flow  conditions, 
where  either  boundary  layer  was  already  turbulent,  should  be  avoided.  With 
the  hot-film  anemometer  held  within  the  boundary  layer  just  upstream  of 
the  T.E.,  the  film  output  was  observed  as  the  stagnation  pressure  was 
changed.  This  qualitative  measurement  is  frequently  used  as  a  rapid  means 
of  ascertaining  the  presence  of  turbulence  (Reference  17).  It  is  based  on 
the  ability  of  this  sensor  to  produce  a  largt-  a.c.  signal  if  immersed  in 
turbulence,  this  signal  decreasing  to  zero  when  the  flow  is  laminar. 

Figure  25  shows  two  such  sets  of  film  signal  vs.  Po  data  for  each 
nozzle,  one  at  the  T.E.  as  indicated,  and  another  at  a  corresponding  point 
on  the  opposite  wall  directly  across  from  the  T.E.  For  the  M  =  2.25  nozzle 
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one  sees  the  signal  begin  from  a  high  level  at  Po2  =  500  mm  Hg,  decreasing 
to  a  low  level  barely  visible  above  the  electronic  noise  when  Po2  =  100  mm. 
The  interpretation  of  these  traces  is  that  for  the  slow-side  nozzle 
turbulence  first  appears  at  about  200  at  the  T.E.,  but  is  not  fully 
established  until  beyond  500  mm.  This  behavior  is  matched  by  the  boundary 
layer  on  the  "ceiling"  of  the  M  =  2.25  nozzle,  as  it  should. 

Without  changing  its  sensitivity,  the  probe  was  next  moved  into 
corresponding  points  in  the  boundary  layer  of  the  Mach  3  nozzle,  where 
the  remaining  two  traces  of  Figure  25  were  obtained.  For  this  case  it 
appears  that  transition  barely  begins  at  Poi  =  600  mm,  while  below  that 
level  the  Mach  3  boundary  layer  is  laminar.  The  conclusion  of  this 
measurement  was  that  Po2  should  he  kept  below  200  mm  Hg  for  all  subsequent 
tests,  while  Po2  could  be  unrestrictedly  set  at  any  level. 

Another  cause  of  constraint  for  Po5  and  Po2  were  the  desire  to  keep 
the  FSL  parallel  to  the  tunnel  axis.  Considering  the  actual  exit  M  for 
the  fast  nozzle  (M2  =  2.8  versus  the  nominal  3.0)  the  two  exit  pressures 
would  be  equalized  if  po2  =  2.3po2. 

The  result  of  these  preliminary  tests  was  the  "design"  of  the  Mode  III 
measurements  at  po2  =  170  mm  and  po2  =  425  mm.  In  addition  to  steadier 
operation  achieved  in  the  tunnel  by  these  two  settings.  Figure  23  shows 
the  possibility  of  observing  transition  some  distance  downstream  of  the 
T.E.  as  desired.  Specifically,  Figure  23  says  that  at  the  T.E.  the  laminar 
FSL  thickness  would  be  0.375  cm,  while  the  expected  transition  thickness 
would  be  about  0.47  cm.  That  is,  the  FSL  was  thus  required  to  first  thicken 
before  it  became  turbulent,  and  transition  was  thus  not  expected  for  some 
distance  behind  the  T.E. 

6.3.2  Mode  III  Nozzle  Flow 


All  subsequent  tests  described  in  this  report  were  done  with  Mode  III 
at  poi  =  425  mm  Hg,  po2  =  170  mm  Hg  and  To2  =  To2  =  605°  R.  The  quality 
of  the  flow  emerging  from  the  two  nozzles  can  be  seen  on  Figure  26.  The 
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uniformity  of  the  velocity  at  the  exit  plane,  which  is  better  than  Vi, 
attests  to  the  care  in  fabricating  the  nozzle  contour. 

Figure  27  shows  the  static  pressure  distribution  in  the  flow  direction 
along  the  ceiling  of  the  "slow"  nozzle,  along  the  floor  of  the  fast 
nozzle,  and  also  along  the  FSL  between  the  two  streams.  The  only  note¬ 
worthy  feature  here  is  an  apparent  longitudinal  pressure  gradient  along 
the  FSL,  amounting  to  (dp/dx)/p  -  0.04  cm’1  (in  terms  of  the  FSL  thickness 
h,  (h/p)  (dp/dx)  =  0.01)  due  to  slight  nozzle  underexpansion.  As  we  shall 
see  later,  however,  flow  quantities  combined  to  produce  much  smaller 
gradients  in  the  other  FSL  edge  properties. 

The  boundary  layer  thicknesses  at  the  T.E.  are  shown  on  Figure  28. 

At  po i  =  425  mm  Hg,  6i  =  0.181  cm,  and  at  po2  =  170  mm  Hg,  S2  =  0.216. 
Although  these  do  not  compare  individually  any  better  than  expected  with 
Low's  theory  (also  shown  on  Figure  28)  their  sum  of  0.397  cm  is  in  almost 
perfect  agreement  with  the  expectation  of  0.396  cm  ->wn  on  the  predictions 
of  Figure  23. 

The  momentum  thickness  at  the  T.E.  is  shown  on  Figure  29.  Some 
reservations  about  the  measured  values  of  02  =  0.0128  cm  and  02  =  0.0158 
are  expressed  at  this  time  because  of  resolution  problems  with  the  pitot 
probe. 

6.3.3  Edge  Properties  and  FSL  Thickness 

Figures  30  through  34  show  the  development  of  the  edge  properties  on 
either  side  of  a  3-inch  (7.5  cm)  length  of  FSL,  including  the  thickness  h, 
the  Mi  and  M2  and  also  p,  u  and  Re^.  The  noteworthy  feature  of  these 
findings,  which  are  needed  to  pin-point  the  transition  Reynolds  number,  is 
that  although  p  has  a  clear  gradient,  the  other  properties  vary  only 
negligibly  along  the  FSL  edge,  as  can  be  seen  from  the  following  Table. 
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TABLE  IV.  Variation  of  Edge  Properties 


p  (mm) 

u  (m/sec) 

M 

Ref  (cm' 

Fast  Side  (1 ) 

Minimum 

13.1 

651 

2.88 

37,400 

Maximum 

16.4 

657 

2.93 

40,000 

Max.  Deviation 

3.3 

6 

0.05 

2,600 

Average 

14.7 

654 

2.9 

38,700 

Slow  Side  (2) 

Minimum 

-- 

577 

2.25 

20  ,400 

Maximum 

-- 

589 

2.32 

21 ,700 

Max.  Deviation 

-- 

12 

0.07 

1  ,300 

Average 

-- 

583 

2.28 

21 ,100 

According  to  this  Table  the  maximum  streamwise  variations  in  u,  M  and 
Re',  are  21,  3%  and  7%  respectively.  Also  note  that,  the  A  as  defined  in 
eq.  (5)  is 

A  =  ^-T-u-2  =  0.0584 
ui  +  u2 


6.3.4  Flow  Profiles  Across  the  Shear  layer 

Figures  35  and  36  show  typical  profiles  across  the  FSL.  Of  these  the 
most  noteworthy  are  the  static  pressure  and  total  temperature  profiles  of 
Figure  35.  The  former  demonstrates  the  uniformity  of  pressure  across  the 
flow,  while  the  latter  shows  that  the  slower  stream  is  about  10°  R  cooler 
than  the  faster  stream.  This  can  be  attribted  to  the  design  of  the 
secondary  channel  through  which  small  heat  losses  were  possible.  Note  that, 
in  anyway,  this  10/600  =  2%  fractional  To  difference  has  only  negligible 
differences  in  the  transition  prediction. 
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The  data  of  Figure  34  indicate  that  nine  sets  of  profiles  were 
taken  along  the  FSL  in  the  region  0<x<3"  from  the  T.E.  As  will  be  shown 
below,  these  were  taken  while  the  FSL  was  still  laminar.  In  the  next 
section,  use  of  these  data  will  be  made  in  fixing  the  FSL  transition 
Reynolds  number  based  on  thickness.  The  data  of  Figures  35  and  36  are 
also  useful  for  future  comparison  with  the  laminar  FSL  theory  described  in 
the  Appendix.  The  velocity  profiles  of  Figure  36  are  especially  useful 
in  confirming  that  the  "wake"  component  is  still  present  for  a  distance  of 
3"  (7.5  cm),  which  corresponds  to  a  distance  of  about  260  total  momentum 
thicknesses  at  the  T.E.  This  is  an  important  finding  indicating  that  the 
flow  profile  has  not  yet  achieved  the  shape  on  which  the  FSL  transition 
theory  is  based,  so  that  agreement  with  the  latter  is  not  necessarily 
expected. 

The  mean  flow  field  in  the  shear  layer  is  discussed  in  greater  detail 
in  Reference  18. 

6.3.5  Transition  to  Turbulence 


The  evidence  as  regards  transition,  which  was  provided  by  the  Schlieren 
system,  is  shown  on  Figures  37,  38,  and  39.  The  former  two  show  the  FSL  as 
the  ratio  P  /P  remains  at  the  "design"  level  of  425/170  =  2.5  but  P 
itself  is  varied  over  the  available  range  200  <  P^  <  600  mm  Hg.  According 
to  the  thinking  expressed  in  Secion  2  of  this  report  transition  is  Reynolds- 
number-controlled,  so  that  for  fixed  X  the  transition  zone  should  move  along 
the  FSL  as  P  is  varied. 

Ol 

This  is  the  case,  according  to  Figures  37  and  38.  The  former  is  a 
series  of  Schlieren  photos  with  continuous  exposure,  while  the  latter  is  a 
a  series  of  spark  photos.  The  onset  of  transition  in  continuous  exposure 
appears,  as  always,  as  a  widening  of  the  FSL  and  a  loss  of  "black-white 
contrast"  due  to  the  decrease  of  density  gradients.  In  the  spark  photos 
a  considerable  amount  of  structure  is  visible  in  the  turbulent  FSL.  In 
both  types  of  photos  the  transition  location  is  quite  clear  and  it  moves 
downstream  as  P  is  decreased. 

Ol 
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Figure  39  shows  spark  and  continuous-exposure  photos  for  the  specific 
conditions  of  the  Mode  III  results  discussed  in  the  previous  Section,  i.e., 

P ^  =  425  mm  Hg  and  Pq  =  175  mm  Hg.  For  these  conditions  it  appears  that 

transition  occurs  at  a  distance  of  2.5"  ±  0.2"  (6.35  ±  0.5  cm)  from  the 
trailing  edge. 

Transition  evidence  in  much  more  detailed  form  was  also  collected 
with  the  hot-film  anemometer,  with  results  shown  on  Figures  40  through  42. 
This  study  exploits  certain  tell-tale  characteristics  of  the  hot-film 
output  behavior  in  geometric  and  Fourier  space.  In  this  instance  a  "map" 
of  the  film  output  in  space  was  first  made  up  by  traversing  the  hot  film 
normal  to  the  FSL  at  several  positions  x  along  a  4-inch  (10  cm)  length,  at 
regular  intervals  away  from  the  T.E.  During  each  traverse  the  output  was 
continuously  recorded  on  seven-track,  1/2-inch  analog  magnetic  tape  at 
120  i.p.s.,  with  a  frequency  response  of  400  khz.  This  method  allows 
careful  review  and  analysis  by  playing  back  the  tape,  after  the  end  of 
the  tunnel  run,  at  reduced  speed. 

Plotted  on  Figure  40  is  the  r.m.s.  wideband  (f  <  400  khz)  hot-film 
output  vs  x  and  lateral  distance  y,  in  the  range  0  <  x  <  10  cm  (or  0  <  x  <  4" 
Considering  the  electronic  noise  magnitude  shown  and  the  origin  ("zero") 
for  each  trace,  it  is  seen  that  a.c.  signals  are  not  obtained  from  the 
FSL  until  x  =  1.5"  from  the  T.E.  Beyond  x  =  1.5"  the  wideband  signal 
becomes  evident  near  the  FSL  center,  indicating  some  sort  of  oscillating  flow 
The  signal  begins  decreasing  toward  x  =  4". 

A  second  "pass"  over  the  same  data  was  next  made  by  looking  at  the 
output  spectrum  each  time  the  film  found  itself  at  the  position  giving  the 
maximum  intensity  wideband  output  (the  "peaks"  of  Figure  40).  Twelve  such 
spectra  are  shown  on  Figure  41.  Each  spectrum  shown  stretches  over  the 
frequency  range  0  <  f  <  320  khz.  Again,  no  oscillations  are  observed  until 
x  =  1.5",  beyond  which  a  strong  oscillation  centered  at  93  khz  is  seen  to 
develop.  The  oscillation  increases  until  it  reaches  a  maximum  at  about  3.25" 
beyond  which  it  begins  decaying. 
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The  93-khz  signal  thus  represents  an  instability  growing  to  detectable 
proportions  by  x  =  1.5"  and  reaching  a  maximum  at  x  =  3.25".  In  fact, 
careful  scrutinizing  of  the  Schlieren  pictures  (Figures  38  and  39)  confirm 
that  such  a  wave- like  phenomemon  is  present  in  the  FSL.  The  spectra  of 
Figure  41  provide  additional  information  on  the  instability.  For  examp  e, 
if  we  assume  that  the  waves  convect  with  the  stream  speed  (u  =  620  m/sec. 
according  to  Table  IV)  then  their  length,  both  absolute  and  relative  to 
the  thickness  h  (cf.  Figure  30)  can  be  computed: 

wavelength  =~  =  0.67  cm  =  1.6  h 

Besides  the  inherent  instability  shown  on  Figure  41,  the  high-frequency 
content  (f>100  khz)  of  the  spectra  also  begins  increasing  at  large  x.  As  already 
discussed  in  connection  with  the  Mode  I,  II  results  (Section  6.2),  this 
represents  the  actual  turbulence  becoming  active.  To  study  this  growth 
separately,  the  traverses  shown  on  Figure  40  were  repeated  once  more, 
except  that  now  only  the  signal  at  f  =  175  khz  (thought  to  be  representative 
of  the  turbulence  content)  was  plotted  against  y.  The  results  are  shown 
on  Figure  42.  The  increase  of  these  high-frequency  signals  as  x  increases 
is  obvious. 

Figure  43  compares  the  downstream  history  of  the  wideband  oscillations 
as  well  as  those  at  93  and  175  khz.  It  is  seen  that  for  the  first  2-2.5  cm 
beyond  the  T.E.  the  FSL  is  completely  quiescent  (this  represents  a  length 
of  about  5-7  thicknesses).  In  the  range  3  -  8  cm  all  Fourier  components 
increase  in  magnitude.  The  instability  peaks  at  about  8  cm,  while  the 
turbulence  itself  peaks  about  1  cm  later,  at  9  cm.  The  simultaneous  growth 
of  discrete  oscillations  along  with  the  growth  of  turbulence,  and  the  decay 
of  the  oscillations  before  the  decay  of  the  turbulence  itself  are  common  to 
transition  in  shear  flows. 

Clearly,  and  as  fully  expected,  the  transition  cannot  be  thus  located 
at  a  "point"  x.  From  Figure  44  it  seems  that  the  process  begins  at  x  =  3 
cm  and  ends  at  9  cm,  for  a  factor  of  3  "spread"  in  wetted  length  x.  This 
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range  corresponds  to  a  wetted- length  Reynolds  number 


from  116,000  to  350,000.  As  will  be  seen  in  the  next  section,  however, 
the  thickness  h  varies  so  little  in  3  <  x  <  9  cm  (see  Figure  30),  that  a 
"point"  in  the  Re^y  =  uih/vx  can  indeed  be  obtained. 

Finally,  a  comparison  of  the  Schlieren  with  the  hot  film  observations 
shows  that  the  Xy  determination  by  the  former  (xy  =  2.5"  =  6.3  cm)  lies  in 
the  middle  of  the  transition  zone  determined  with  the  latter.  Partly 
because  it  thus  acts  as  an  average  to  the  range,  the  distance  Xy  =  2.5"  will 
be  henceforth  considered  as  the  measured  transition  "point." 
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Discussion  of  the  Transition  Observations 


The  Mode  III  results  give  the  main  quantitative  indications  of 
transition  location,  which  can  be  discussed  in  the  context  of  the 
asymptotic  theories  presented  in  Section  2.  To  this  end,  the  transition 
Reynolds  number  Re^  was  formed  from  the  data,  using  the  average  fast-side 
unit  Reynolds  number  Ref  =  38,700  cm'1,  and  the  indicated  dimensionless 
speed  ratio  \  =  0.0594  (see  Section  6.2.3).  The  FSL  thickness  h  where 
transition  occurs  is  obtained  from  Figure  43.  Using  h  =  0.415  cm  at 
x  =  2.5",  one  obtains  Re^  =  15,500  by  which  the  datum  on  Figure  44  is 
plotted. 

Figure  44  also  plots  the  FSL  transition  expectation  as  obtained  from 
eq.  (4)  and  Figure  2,  for  the  average  observed  Mi  =  2.9  (see  Table  IV). 

The  asymptotic  limit  (X  =  o)  for  a  wake-type  flow  is  also  plotted,  first 
by  computing  the  requisite  Rew: 

Re  =  2Re(  (0!  +  02)  =  38,700  x  2  x  0.0286  =  2200 

W 

and  then  finding  Re^y  »  22,500  from  eq.  (4)  or  Figure  2. 

Several  critical  comments  can  be  made  about  the  results  shown  on 
Figure  44.  First,  it  appears  that  the  agreement  between  data  and  predictions 
is  very  good,  especially  considering  the  nature  of  the  theory  and  the 
difficulty  in  defining  transition  onset.  Secondly,  it  can  be  said  that,  if 
the  theoretical  curves  are  taken  literally,  the  25%  difference  between  the 
test  datum  and  the  prediction  represents  the  expected  departure  from  the 
theory  due  to  the  "hybrid"  flow  profile  which  is  neither  a  wake  nor  a  FSL 
(see  profiles  at  x  =  2"  and  3"  in  Figure  36). 

Third  it  can  be  speculated,  as  already  done  ealier  on  Figure  6,  that 
the  actual  transition  Reynolds  number  curve  joins  smoothly  the  >  =  o  and 
X  =  1  limit  via  the  dotted  portion  shown  as  "speculative"  on  Figure  44.  If 
this  is  true,  then  the  datum  obtained  bridges  an  important  gap  in  the  under¬ 
standing  of  transition  as  viewed  by  the  present  theory  of  Section  2. 
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Finally,  it  must  be  kept  in  mind  that  h  grows  slowly  with  x,  which  is 
true  for  almost  all  boundary-layer  flows.  If  a  prediction  that  something 
will  occur  at  a  certain  value  of  h  carries  along  a  probable  error,  this 
error  is  greatly  magnified  in  terms  of  x.  Thus,  even  if  the  present 
theory  is  capable  of  good  transition  estimates  in  terms  of  FSL  thickness, 
it  can  predict  the  distance  of  transition  from  the  T.E.  only  to  the  extent 
to  which  reliable  predictions  for  h  (x)  exist. 

Quite  separate  from  the  transition-location  issue  is  that  of  the  role 
of  instabilities  in  the  shear  layer.  Although  instabilities  are  a  natural 
precedent  of  transition  in  free  shear  flows,  the  attention  they  have 
recently  attracted  is  more  due  to  their  potential  as  a  cause  of  large-scale 
mixing.  If  one  was  to  judge  from  the  evidence  at  low  speeds  alone,  such 
"laminar  mixing"  could  be  of  great  importance,  say,  to  laser  devices.  The 
present  experiments,  however,  indicate  conditions  for  which  unstable  vortex 
or  wave  structures  develop  only  after  some  distance  downstream  of  the  T.E.; 
and  once  they  do,  their  amplitude  is  still  no  greater  than  the  laminar  FSL 
thickness  itself.  It  is  therefore  submitted  that  for  such  cases  large 
periodic  motions  of  the  FSL  do  not  seriously  contend  to  the  mixing  process. 
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Conclusions 


The  following  was  learned  from  the  work  performed  so  far: 

1)  The  threshold-based  transition  theory  for  free  shear  layers  was 
found  to  need  modification  when  the  shear  layer  is  formed  by  a 
partition  initially  separating  the  two  streams.  By  recognizing 
the  existence  of  the  boundary  layers  shed  from  the  partition, 
finite  transition  Reynolds  numbers  were  obtained  for  the  entire 
range  of  speed  ratios  when  such  modifications  were  made. 

2)  A  closed-form  analytic  solution  for  the  FSL  flow  from  the  T.E. 
onward  was  constructed,  which  appears  to  be  as  accurate  as  exact 
(numerical)  solutions,  and  which  is  much  easier  to  use  than  the 
latter. 

3)  Shear-layer  transition  data  obtained  for  M;  =  3,  1.6  <  M  <  2.2  and 
0.05  <  X  <  0.15  were  completely  consistent  with  the  theory  mentioned 
under  (1),  above.  Specifically, 

(a)  When  the  theory  predicted  turbulence  for  a  FSL  shed  off 
the  T.E.  while  still  laminar,  the  layer  became  turbulent 
immediately  past  the  T.E. 

(b)  When  the  theory  predicted  that  the  FSL  would  become 
turbulent  only  after  some  distance  past  the  T.E.  elapsed, 
this  was  found  to  be  true  experimentally.  For  X  =  0.058 
the  observations  gave  a  transition  Reynolds  number  of  16,000 
which  is  within  25%  of  the  prediction. 

4)  Laminar  shear  layers  were  obtained  which  were  free  of  periodic 
oscillations.  Such  oscillations  were  in  that  case  only  present  as 
pre-transitional  instabilities,  and  were  of  amplitude  and  duration 
small  enough  to  make  them  incapable  of  laminar-laminar  mixing  by 
themselves. 
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Figure  7.  Transition  on  FSL's  shed 
from  partitions  at  Mi  =  3. 
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Figure  11.  SWT  modification  for 

parallel-stream  mixing. 


configuration.  Note 
probe  installation. 


SIS 


16.  Boundary-layer  growth  in 
the  SWT  Mach  3  nozzle. 
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Ill  flow.  Spark  photo  conditions  are 
(middle)  and  286  (bottom)  mm  Hq. 


FLOW  ACROSS  NOZZLE  EXITS 


both  nozzles. 


DISTANCE  FROM  M=3  THROAT  (INCH) 


Figure  27.  Pressure  distribution 
in  the  nozzle  exit 
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Figure  34.  Center-plane  velocity 
and  temperature 
changes.  Mode  III. 
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continuous  (top)  and  spark 
(bottom)  Schlieren 
photography . 


Figure  42.  The  hot-film  output  at 

f  =  175  khz  across  the  FSL 
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APPENDIX  A 


THE  LAMINAR  SHEAR  LAYER  STARTING 
AT  THE  TRAILING  EDGE 


1 .  Summary 

A  theoretical  solution  is  presented  of  the  two-dimensional  laminar 
shear  layer  beginning  at  the  trailing  edge  of  a  partition  separating 
two  homogeneous,  steady  but  dissimilar  parallel  flows.  Dissimilarity 
applies  not  only  to  the  stream  velocities  and  temperatures  but  also  to 
the  initial  boundary  layer  thicknesses.  The  validity  region  of  the 
solution  extends  from,  to  infinite  distances  downstream  of,  the  trail¬ 
ing  edge. 

2.  Introduction 


The  gas-dynamic  and  chemical  laser  technology  has  motivated  a  fresh 
look  at  laminar  free  shear  layers  (FSL)  shed  of  the  trailing  edge  (T.E.) 
of  partitions  separating  parallel  flows.  As  posed,  the  problem  usually 
addresses  supersonic  flows,  usually  in  the  range  2  <  M  <  6.  Furthermore, 
the  technology  requirements  focus  on  characteristic  wetted  lengths  which 
are  too  small  for  turbulent  flows;  the  surface  boundary  layers  are  there¬ 
fore  most  likely  laminar  as  is  the  resulting  initial  FSL.  Figure  A.l 
sketches  the  flow  profiles  one  intuitively  expects  to  see  downstream  of 
the  partition  T.E.  For  simplicity  the  partition  is  assumed  infinitely 
thin,  and  the  flows  homogeneous  (chemically  identical  and  inert),  steady 
and  laminar. 

Subscripts  1  and  2  on  Figure  A.l  refer  to  the  fast  and  slow  stream  sides 
(ui>u2)  respectively.  The  flow  properties  u,  M,  Tq  etc.  are  allowed  to 
differ  between  the  two  sides.  Differences  are  also  allowed  in  the  initial 
boundary  layer  thicknesses  S  and  0.  The  problem  posed  is  the  prediction 
of  the  profiles  of  u,  T,  o  etc.  as  a  function  of  x*  and  y*  (the  physical 
longitudinal  and  lateral  distances)  or  their  non-dimensional  counterparts. 
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The  cusp-like  initial  flow  profile  at  the  T.E.,  predictable  from 
laminar  boundary  layer  theory,  will  eventually  become  the  familiar  FSL 
profile  far  from  the  trailing  edge.  This  asymptotic  behavior  is  due 
to  the  inevitable  diffusion  of  the  "wake  component"  caused  by  the  boundary 
layer  flow  on  the  partition.  Because  of  its  invariability  with  x.  the 
final  profile  is  here  termed  the  "equilibrium"  profile;  the  term  "non¬ 
equilibrium"  is  used  to  denote  the  zone  between  the  T.E.  and  the  equilib¬ 
rium  zone,  where  the  wake  component  is  visible  and/or  dominant.  The 
laminar  non-equilibrium  zone  is  addressed  here  because  of  its  potential 
importance  to  mixing,  stability  and  transition. 

While  the  equi 1 ibrium-zone  profiles  have  attracted  much  successful 

theoretical  study,  there  is  a  surprising  shortage  of  analytical  attention 

to  the  non-equilibrium  flow.  This  is  perhaps  due  to  the  conviction  that 

solutions  can  be  obtained  numerically.  Much  more  useful  however  would  be 

an  analytic  solution  in  closed  form,  without  restrictions  on  the  range  of 

x,  and  capable  of  handling  the  full  range  of  compressibility  (0  <-M  <  =°), 

asymmetry  (0^0,,  ui/uO  and  heat  transfer  (T  fl  )  conditions.  Such  a 

0 1  0  2 

solution  is  presented  below. 

Method  of  Analysis 

The  task  involved  is  to  obtain  a  closed-form,  analytic  expression  for 
the  entire  flow  field  of  Figure  A.l,  complete  with  compressibility,  heat 
transfer  (T  j*T  )  an<1  arbitrary  values  of 


To  offset  the  magnitude  of  the  task,  the  pressure  is  assumed  everywhere 
constant.  Also,  it  is  clear  that  the  precise  shape  of  the  initial  boundary 
layer  profile  will  be  insignificant  to  the  downstream  flow. 

Two  earlier  approaches  to  this  problem  deserve  mention,  both  addressing 
"base  flow"  problem  (r=0)  in  vogue  in  the  1960's.  Dennison  and  Baum 
(keference  A.l)  found  a  numerical  solution,  which  is  presumably  exact. 

Kubota  and  Dewey  (Reference  A. 2)  approached  the  same  problem  with  the 
momentum-integral  method;  the  main  disadvantage  of  this  method  (also  used 
subsequently  by  Ortwerth  and  Moeny  (Reference  A. 3)  as  a  solution  for  cases 
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with  r  f  0)  is  the  existence  of  discontinuities  in  the  flow  profiles. 

Both  the  Dennison-Baum  end  Kubota-Dewey  solutions  demonstrated  agreement 
with  the  large-x  FSL  profile  of  Chapman  (Reference  A. 4).  The  Dennison- 
Baum  and  Chapman  solutions  will  be  of  further  use  here,  the  latter  as  a 
direct  input  to  the  present  theory,  and  the  former  as  a  gage  of  its  validity. 

The  present  scheme  is  based  on  the  possibility  of  superposing  linear 
flow  fields,  and  in  broad  outline  follows  this  procedure:  it  is  recognized 
that  the  flow  at  hand  consists  of  two  separate  components,  the  wake  flow 
generated  by  the  two  merging  boundary  layers  on  one  hand,  and  the  usual 
FSL  profile  joining  the  two  parallel  streams  on  the  other.  The  shear- 
layer  component  can  be  visualized  as  starting  at  the  T.E.  and  it,  of  course, 
dominates  the  flow  at  far  downstream  distances.  The  wake  component  is  most 
visible  at  and  shortly  downstream  of  the  T.E.,  but  it  has  little  presence 
far  downstream  because  of  lateral  diffusion.  The  analysis  then  consists  of 
the  linear  addition  of  the  two  flow  fields,  that  is  of  following  the  wake 
component  as  it  develops  on  a  "baseline"  contorted  to  conform  to  the  shear- 
layer  component. 

The  idea  promoting  the  above  method  comes,  essentially,  from  experi¬ 
mental  observations  dealing  with  wakes  growing  in  non-uniform  flows.  In 
low-speed  flows  Eskinazi  (Reference  A. 5)  showed  that  the  turbulent  wake  of 
a  cylinder  growing  across  a  pipe-flow  profile  was  not  affected  by  lateral 
velocity  gradients.  Kubota  and  McCarthy,  using  the  laminar  wake  of  a  cyl¬ 
inder  in  hypersonic  flow,  actually  demonstrated  that  laminar  wake  similarity 
was  found  when  the  laterally  non-uniform  flow-field,  generated  by  the  bow 
shock,  was  "subtracted  out"  (Reference  A. 6).  Motivated  by  this  finding, 
the  present  author  (Reference  A. 7)  showed  that  when  the  shear-layer  compo¬ 
nent  due  to  the  supersonic  inviscid  flow  generated  by  a  plate  at  incidence 
was  similarly  "subtracted",  the  resulting  wake  profile  conformed  to  similar¬ 
ity  predictions. 

Consider,  therefore,  two  streams  labeled  1  and  2  and  pictured  on 
Figure  A.l  such  that  the  "speed  ratio"  r  =  u2/ui,  flowing  along  x*,  where 
x*  =  0  at  the  T.E.  (starred  quantities  are  physical)  and  use  the  notation 


y'  to  signify  the  compressible-transformed  coordinate  y  divided  by  the 
total  momentum  thickness: 


' 

~  '  e{ +e-t_ 


H* 

(’ ?  I  * 
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Then  the  wake  component  of  the  velocity: 


U.  3  I  — 


U4 


X  = 


f?e0 


Vk 


(A.  1 ) 

(A. 2) 
(A. 3) 


where  "e"  is  used  to  signify  the  "external"  or  "baseline"  stream.  Also 
consider  the  shear  layer  between  the  two  streams, 

Ms  s  rr,  (A-4) 

so  that  the  combined  velocity  field  will  contain  the  characteristic  asym¬ 
metry  parameters  r  and  P  of  its  two  individual  components: 

U  =  U(u  ,  3  U  (  %'j  y'  j  PJ  r)  (A. 5) 

The  two  components  will  be  briefly  discussed  below. 

4.  The  Wake  Component 


The  wake  velocity  field  for  arbitrary  P  has  already  been  provided  by 
this  author  (References  A. 8  and  A. 9).  The  solution  is  analytic,  straight¬ 
forward  and  simple: 

,  r^(P+o\'+  (P+Oy' 

1  !. 


u 


U.  =  1—77  *  —  !  e  (  I  -  E *1  £  ((  P-m)  ^X  '  (A. 6) 
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+  L  )) 

z\T x'  1 


+  e 


.  P+1 /  P+' 

C  O  )  x  -  .p 


This  solution  is  based  on  assumed  exponential  boundary  layer  profiles  and 
utilizes  an  Oseen-type  linearization  (Reference  A. 10).  However,  the  solution 
obeys  all  the  required  initial  and  boundary  conditions,  and  agrees  very 
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well  with  Goldstein's  implicit  wake  solutions  for  P  =  1  (see  Reference  A. 11). 
Furthermore,  the  initial  profile  should  be  unimportant  past  the  T.E.  so 
long  as  the  normalizing  length  is  the  total  boundary- layer  thickness 
(Reference  A. 8). 


5.  The  Shear-Layer  Component 


5 . 1  Chapman's  Solution 


To  complete  the  solution,  the  quantity  u$  appearing  in  eq. 
needed.  One  possibility  is  the  solution  first  given  by  Chapman 
A. 4)  for  r  =  0,  in  the  form  u/ui  vs  the  similarity  variable 


'h 


(A. 4)  is 
(Reference 

(A. 7) 


Since  u/ui(n")  was  computed  by  Chapman  by  numerical  evaluation  of  certain 
integrals,  it  was  found  convenient  here  to  approximate  it  by  a  least- 
squares  curve-fit  polynomial: 


3 


vi 


(A. 8) 


The  variable 
by 


is  related  to  the  x',  y'  of  the  wake  solution  (eq. 


<ffi  *e. 


x/ 


(A. 6)) 
(A. 9) 


where  the  identity  between  y*  and  its  compressible-transformed  y*  is  assumed 
since  eq.  (A. 8)  was  obtained  from  Chapman’s  incompressible  (M=0)  result. 
Furthermore,  it  was  more  convenient  to  replace  n'  by 

(A. 10) 


'  2_'Tx"/ 


=  "1 


so  that  eq.  (A. 8)  was  changed  to 

4 


U. 

u 


(A. 11) 


n  ro 


with  C0  =  0.587,  Ci  =  0.381292,  C2  =  5.7138  x  10'4,  C3  =  -5.238816  x  10~2 , 
C„  =  -5.244608  x  10-J,  C5  *  3.465856  x  10"’,  C5  =  6.070848  x  10'4.  Note 
that  eq.  (A. 11)  is  valid  only  in  -4  <  n  <  2.5;  if  n  >  2.5  then  u  =  ui,  and 
if  n  <  -4  then  u  =  0. 
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To  obtain  the  ug  quantity  from  eq.  (A. 11),  for  input  into  eq.  (A. 4), 
it  was  necessary  to  adapt  Chapman's  result  to  the  case  u2  t  0,  i.e.  for 
0  <  r  <  1.  Since  the  only  relevant  FSL  motion  is  the  relative  one  between 
the  two  streams,  then  Chapman's  profile  of  eq.  (A. 11),  originally  found 
for  r  =  0,  will  give,  when  r  f  0: 


U&  3  =  C!-r)  ^  Cn>?  + 

*  ^  o  L 

which  is  the  desired  result  for  insertion  into  (A. 4). 
appears  on  Figure  A. 2. 


(A. 12) 

The  form  of  eq.  (A. 11 ) 


5.2  Mills'  Solution 


A  second  possibility  for  u$  is  the  well-known  error-function  formula 
first  given  by  Goertler  (Reference  A. 12): 


u 


s 


/  +  r 

T 


i  ~r 


2. 


(A. 13) 


where  o  is  the  appropriate  similarity  variable  combining  the  longitudinal 
and  lateral  coordinates.  This  solution  has  undergone  considerable  scrutiny 
since  its  inception.  More  recently,  a  series  of  papers  by  Lock  (Reference 
A. 12),  Ting  (Reference  A. 13),  Crane  (Reference  A. 14)  and  Mills  (Reference 
A. 15)  has  re-examined  the  problem  by  seeking  the  elusive  "third  boundary 
condition"  to  the  equations  of  motion.  Mills'  paper,  being  the  last  in  the 
series  chronologically,  reviews  the  issue  thoroughly  and  arrives  at  a 
numerical  method  of  computing  the  laminar  shear  layer  for  any  r,  Mj  or  Tq  /T 
The  gist  of  Mills'  paper  is  that  although  eq.  (A. 13)  is  a  good  approximation 
to  the  shape  of  the  profile,  the  value  of  u/ui  needs  adjustment  by  some 
function  of  r  such  that  the  argument  of  the  error  function  is  n  +  a(r),  where 
a(r)  is  a  function  he  tabulates.  His  tabulation  was  used  to  find  the  follow¬ 
ing  curve-fit: 

0.tr)  =?0.lL>H-7t>7  -  0.  0  I  k  'Z'j  S  7r  *0  it  if  r%  q.  3  3  37  f 7  r  3  (A  14) 


Mills  Furthermore  expresses  n  as 

<j  *  — .  t 

'"I  »  (RO  *- 
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Re'  S  ) 


(A. 15) 
(A. 16) 
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where  Re'  is  the  unit  Reynolds  number.  Equations  (A. 14),  (A. 15)  and  (A. 16) 
comprise  the  alternative  shear-layer  velocity  profile  to  be  used  in  (A. 4): 


Us  -  E  ( "n  +•  afrb)  (A.  17) 

2  2.  1 

Equation  (A. 17)  is  compared  with  Chapman's  method  on  Figure  A. 2. 

6.  Combination  of  the  Two  Flow  Components 

In  combining  the  two  fields,  we  assume  consistently  with  the  experi¬ 
mental  findings  quoted  above  that 

(a)  The  wake  component  develops  independently  of  the  shear  component. 

Equation  (A. 6)  presents  no  obstacle  to  this  proposition  if  the  free-stream 
velocity  ug  in  the  definition  of  u  is  taken  to  be  the  shear  velocity  u', 
except  that  the  unit  Reynolds  number  Re'  entering  the  non-dimensional 
coordinate  x'  has  to  be  determined.  As  a  concession  to  the  difference 
between  Re',  and  Re'2,  we  further  propose  that 

(b)  The  unit  Reynolds  number  controlling  the  diffusion  of  the  wake 
component  of  eq.  (A. 6)  is  the  average  between  Re'i  and  Re'2,  as 
defined  in  (A. 10). 


These  two  propositions  enable  us  to  use  the  wake  component  as  it  stands 
eq.  (A. 6)  with  ^  ^ 


'  -  d 

'd  =  - 

0 


/ 

y  s 


x 


Da  - 


5  S 


9 


(A. 18) 


For  the  shear  component  we  assume  that  the  lateral  coordinate  y  is  replaced 
by  its  compressible-transformed  counterpart  y,  which  is  hardle  controver¬ 
sial;  then  the  coordinate  n  in  (A. 9)  immediately  becomes 


1 


2  \Tx' 


(A. 19) 


For  the  shear  component  u$  we  have  a  choice  between  Chapman's  form 
(eq.  (A. 12))  and  Mills'  form  (eq.  (A. 17)).  Note  that  in  both  cases  the 
independent  variable  n  is  given  by  eq.  (A. 9).  However,  a  point  arises  in 
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the  physical  interpretation  of  n  when  the  wake  and  shear-layer  components 
are  combined,  as  will  be  done  below.  Specifically,  while  proposition  (b), 
above,  applies  to  the  wake  as  well  as  to  Mills'  FSL  component  (eq.  (A. 17)), 
it  is  clearly  not  applicable  to  Chapman's  case  when  Re2  =  0.  For  the  moment 
we  will  nevertheless  also  assume  that  the  Reynolds  number  appearing  in 
(A. 12)  is  the  average  given  by  (A. 10)  if  u2  f  0;  this  makes  x'  uniform 
throughout. 


The  addition  (superposition)  of  the  two  components  is  then  explained 
with  the  aid  of  Figure  A. 3.  The  actual  flow  velocity  u  is  related  to  the 
shear-component  velocity  u'  by 

U  =  u'  -Uw  (A. 20) 


In  ordinary  wakes,  the  u  is  usually  the  difference  between  the  edge  or 

W 

external  velocity  ug  and  the  actual  velocity: 

Uw  *  U{  u  (A. 21) 

where  u  is  already  defined  in  eq.  (A. 6).  The  key  to  the  present  approach 

is  to  set  u '  =  u  ,  i.e.  to  superpose  the  wake  flow  on  a  "free  stream" 

6 

which  is  actually  the  shear-layer  profile.  Inserting  u'  =  ug  into  the 
above  two  formulas. 

Us  =  U  '(  I  -  u  )  (A. 22) 


or 


u 

u 


^  (  (  _  U  )  =  usC!  -  U  ) 


(A. 23) 


with  us  supplied  by  (A. 12)  or  (A. 17)  and  u  by  (A. 6).  Thus, 


‘I  ‘-Jv 

The  thermodynamic  variables  can  then  be  obtained  with  the  aid  of  the 
Crocco  relation: 


(A. 24) 


r  +  V  (A. 25, 


from  which  the  temperature  ratio  is 

T,  -T  Ky'j  r,1>i  ) 


(A. 26) 


Once  the  temperature  ratio  is  known,  the  remaining  thermodynamic  parameters 
can  be  calculated  easily. 
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7.  Illustration  and  Discussion  of  the  Solutions 


7.1  Computation  Method 

Three  simple  BASIC  language  problems  have  been  written  to  compute  the 
flow  according  to  eq.  (A. 23).  Program  SHEAR4  requires  inputs  of  P,  r.  Mi, 
T.  /T  ,  and  x'  and  produces  profiles  of  temperature,  velocity  etc.  as  a 

□  2  0 1  „ 

function  of  y  (the  required  y'  interval  is  an  input).  This  program  uses 
Mills'  shear  Layer  Model,  eq.  (A. 17).  Program  SHEAR5  requires  the  same 
inputs,  and  produces  only  profiles  vs.  y'.  Th.jse  profiles  are  of  the  (a) 
Mills'  FSL  component  from  eq.  (A. 17),  (b)  the  complete  flow  from  eq.  (A. 23) 
using  Mills'  u$,  (c)  the  Chapman  FSL  component  from  eq.  (A. 12),  and  (d)  the 
complete  flow  from  eq.  (A. 23)  using  Chapman's  component.  Finally,  program 
SHEAR6  computes  properties  only  along  the  y'  =  0  plane. 

These  programs  are  shown  on  Tables  I,  II  and  III. 
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8.  Results 


Computations  with  equations  (A. 12),  (A. 17)  and  (A. 23)  have  been  carried 
out  and  are  shown  on  Figures  A. 4  through  A. 10.  Typical  questions  answered 
by  these  computations  are  (a)  how  does  the  present  theories  agree  with  other 
methods,  (b)  how  do  results  using  Mills'  u$  differ  from  those  using  Chapman's 
u$,  (c)  how  is  the  "regular"  FSL  profile  get  distorted,  and  for  how  long, 
by  the  shed  boundary  layers? 

(a)  Comparison  with  other  theories  is  made  difficult  by  the  absence 
of  parallel  methods  addressing  the  issue  as  widely  as  done  here.  A  good 
test  can,  however,  be  made  in  the  case  of  r  =  0,  i.e.  the  "base  flow" 
problem  addressed  by  Kubota  and  Dewey  and  by  Dennison  and  Baum.  Figure  A. 5 
shows  the  quantity  u/ui  (y'  =  0)  for  this  case.  The  agreement  of  the  present 
method  (using  Chapman's  ug)  with  the  Dennison-Baum  (exact)  solution  is  very 
good,  and  certainly  tjtter  than  Kubota  and  Dewey  found  by  the  momentum- 
integral  method.  Note  that  the  present  theory  is  much  easier  to  apply  than 
either  of  these  other  two  methods. 

(b)  Differences  in  the  results  are  found,  throughout  the  attached 
Figures,  between  results  obtained  with  Mills'  u$  and  those  found  by  Chapman's 
u$.  In  general,  the  Chapman  approach  gives  lower  velocities.  The  difference 
tends  to  disappear  when  r  is  close  to  unity,  however  (e.g.  compare  Figure 

A. 8  with  Figure  A. 9). 

(c)  Whether  the  flow  profile  "looks  like  a  shear  layer"  or  a  "wake" 
can  be  directly  seen  from  eq.  (A. 23)  where  the  influence  of  the  wake  compo¬ 
nent  appears  as  the  factor  (1  -  u).  For  example  if  we  require  that  the 
regular  FSL  profile  u$  is  not  distorted  by  more  than,  say,  5 %,  we  find  the 
x'  at  which  u  has  decreased  to  0.05  on  the  axis.  It  turns  out  that  this 
distance  is  very  long,  of  order  x"  =  100.  Maximum  x'  attainable  in  laser 
cavities  is  of  order  1-10,  and  thus  it  follows  the  laminar  FSL's  in  such 
cavities  require  consideration  of  the  wake  component  as  done  herein.  Note 
that  the  influence  of  the  wake  component  cannot  be  always  gaged  by  the 
presence  or  absence  of  the  "trough"  or  "bucket"  in  the  flow  center.  For 
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example.  Figure  A. 10  shows  that  by  x'  =  0.1  the  trough  has  been  "smeared 
out"  by  the  flow  and  yet  the  flow  velocity  is,  in  places,  half  of  what  it 
would  be  without  the  initial  boundary  layers. 
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TABLE  IA.  THE  SHEAR4  PROGRAM 


tSCH  <  <  )-</  (2»SW<» 


TABLE  IA.  (Continued) 
PRINCIPAL  SHEAR4  PROGRAM  GLOSSARY 


Program 

Symbol 

First  Appearing 
on  Line  Number 

Is  the 
Algebraic 

Meaning 

X 

20 

X' 

x*/0'Re^ 

M 

22 

Mi 

-- 

T 

24 

VTo, 

-- 

PI 

30 

P 

9i/02 

R 

40 

r 

Ui/Ui 

B1....B4 

120-150 

- 

Coefficients  of  a(r) 

A 

170 

a(r) 

-  “ 

P 

230 

n 

Constant  in  the 
approximation  to  Erf 

A1....A5 

240-280 

A i , ...  A  s 

Constants  in  the 
approximation  to  Erf 

Y 

310 

y ' 

yVO 

VI 

320 

(p+1  )*r~x'  +  ti_ _ 

X  1 

IfT* 

FI 

400 

Xlexp((P+D2x'  + 

(P+i )y ') 

X2 

410 

y  ■*  P+1  y~~%* 

t7r-~  'r 

Argument  of  2nd  Erf 
appearing  in  u 

F2 

490 

Second  term  in  u 

U 

500 

u 

Wake  component  of 
velocity 

H 

510 

n 

y72i •r*r~ 

K 

590 

Erf (o+a(r)) 

Erf  in  u$ 

U1 

600 

us 

FSL  component  of 
velocity 

U2 

650 

u/Ui 

-- 

T4 

680 

T/T  i 

-- 

M2 

690 

M 

-- 

TABLE  IA.  (Continued) 
SHEAR4  PROGRAM  LANDMARKS 


LINES 

OPERATION 

120-190 

Computation  of  a(r),  eq.  (8. A) 

230-290 

Definition  of  constants  in  approximate  formula 
for  Erf 

310 

Loop  start  (sets  y'  intervals  for  printout) 

320-390 

Computation  of  complementary  Erf  in  the  first 
u  term 

400 

Computation  of  first  u  term 

420-480 

Computation  of  Erf  in  second  u  term 

500 

Final  computation  of  U 

530-590 

Computation  of  the  Erf  in  us 

600-620 

Final  computation  of  u$ 

630-680 

Computation  of  T/Ti 

650 

Computation  of  u/ui 

690 

Computation  of  M 
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TABLE  IA.  (Continued) 


ACTIVATION  COMMANDS  FOR 
SHEAR4,  -5  AND  -6 


LINE 

COMMAND 

20 

Enter  x' 

22 

Enter  M] 

24 

Enter  T  /T 

Ot  Oi 

30 

Enter  P 

40 

Enter  r 

310 

Enter  y'  range  and  intervals 

305* 

Enter  x"  range  and  intervals  by  typing  in  N 
range  and  intervals  where 

*  SHEAR6  only 
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TABLE  IIA.  THE  SHEARS  PROGRAM 
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TABLE  1 1 A .  THE  SHEAR6  PROGRAM 


Figure  1.  Nomenclature  and 

Definitions  of  the  FSL 
Problem 


Figure  4.  FSL  for  r  =  0  for  Mills 

(left)  and  Chapman  (right) 
models. 
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Figure  6.  Typical  FSL  velocity 
profiles  (r  =  0.1 ) 


10  -8  -4  0  4  8  10 

Figure  7.  Typical  FSL  velocity 
profiles  (r  =  0.5) 
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Effect  of  wake  component 
on  FSL  flow  (typical ). 


